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PATENT 



*§4 OBt f0? / IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

x ... . Herbert Confirmation onKQ 

Applicant: BRUDERet al. No.: 9058 

Application No.: 10/667,475 Examiner: John Corbett 

Filing Date: f^ mber23. Group Art Unit: 2882 
Title: 

Attorney Docket: 32860-000622/US 



2003 

METHOD OF CREATING IMAGES IN COMPUTER 
TOMOGRAPHY (CT), AND CT DEVICE 



DECLARATION UNDER 37 C.F.R. g 1.131 



COMMISSIONER OF PATENTS 
P.O. BOX 1450 

ALEXANDRIA, VA 22313-1450 

Sir: 



We, Herbert Bruder, Thomas Flohr, Annabella Rauscher-Scheibe, Karl 
Schwarz and Karl Stierstorfer, declare as Inventors of the subject matter 

claimed in the above-identified application that: 



1. This Declaration is submitted as evidence that the subject 
matter disclosed and claimed in U.S. Patent Application No. 10/667,475 
was invented by us prior to July 10, 2002, which is the effective filing date of 
U.S. Patent NO. 7.245.755 Bl ("Pan") entitled "Algorithm for Image 
Reconstruction and Image Noise Analysis in Computed Tomography." 

2. We conceived, and are the named inventors of, the subject 
matter disclosed and originally claimed in U.S. Patent Application No. 
10/667,475. 



AppUcation No.: 10/667,475 
Attorney Docket No.: 32860-000622/US 

3. U.S. Patent AppUcation No. 10/667,475 is currently assigned to 
Siemens Aktiengesellschaft (hereinafter interchangeably referred to as "the 
company") as evidenced by the assignment recorded on Reel 014973, Frame 
005 1 in the U.S. Patent and Trademark Office. 

4. The company is a company doing business in the Germany, a 
World Trade Organization fWTO) member country. 

5. Prior to July 10, 2002, we conceived a method of creating 
images in computer tomography (CD and a CT device in Germany, which at 
the time of conception was a WTO member country. (A copy of an internal 
Invention Disclosure Form and English translation thereof are attached as 
ExhibitAJ. 

6. Hie Invention Disclosure Form was submitted on April 30, 
2002, handed over to a supervisor, Dr. Stephan Schaller, on May 7, 2002, 
and was received in the Patent Department for the company on May 23. 
2002. 

7. The subject matter disclosed and claimed in U.S. Patent 
AppUcation No. 10/667,475 was disclosed in the Invention Disclosure Form 
entitled "Gewichtungsverfahren fur 3D-Ruckprojeckttonsverfahren fur die 
Mehrzetfer-Spiral-CT shown in Exhibit A. 



Application No.: 10/667,475 
Attorney Docket No.: 32860-000622/US 

8. Prior to July 10, 2002, outside counsel (Dr. Annin Binder) was 
requested to prepare a patent application based on our Invention (A copy of 
the order letter dated Jury 3, 2002 and an English translation thereof are 
attached as Exhibit B). 

9. On August 28, 2002, Dr. Binder forwarded the draft patent 
application to one of us, Dr. Karl Stierstorfer, for review (A copy of a letter 
from Dr. Binder to Dr. Stierstorfer and an English translation thereof are 
attached as Exhibit C). 

10. After reviewing the draft application, comments regarding the 
draft application were sent to Dr. Binder (Copies of emails dated August 29, 
2002 and September 4, 2002 and English translations thereof are attached 
as Exhibits D and E, respectively). 

11. On September 23. 2002, German Patent Application 10244181.2 
(the *"181 Application), which conformed substantially to the draft forwarded 
from Dr. Binder to Dr. Stierstorfer, but with the comments shown in 
Exhibits D and E, was filed with the German Patent Office (A copy of the 
'93i2 Application and an English translation thereof are attached as Exhibit 
F). 

12. On September 23, 2003, the above-identified application was 
filed as U.S. Patent Application No. 10/667,475, and claimed priority under 
35 U.S.C § 119 to the '181 Application, Hie U.S. application is an English 



Application No.: 10/667,475 
Attorney Docket No.: 32860-000622/US 



translation of the filed German Patent Application, slightly modified to place 

i 

it In better form for U.S. practice. 

13. We hereby declare that all statements made herein of our own 
knowledge are true and that all statements made on information and belief 
are believed to be true; and further that these statements were made with 
the knowledge that willful false statements and the like so made are 
punishable by fine and imprisonment, or both, under Title 18, United States 
Code, § 1001 and that such willful false statements may jeopardize the 
validity of the application or any patent issued thereon. 



Herbert BRUDEJ* Dated 



Thomas FLOHR 



Dated 




Annabella RAUSCHER - SCHEIBP 



Dated 



Application No.: 10/667,475 
Attorney Docket No. : 32860-000622/US 

Karl SCHWARZ Dated 



Karl STIERSTORFER Dated 
89 J 495.1 
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Attorney Docket: 32860-000622/US 
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No.: 
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Examiner: John Corbett 

Group Art Unit 2882 



2003 

METHOD OF CREATING IMAGES IN COMPUTER 
TOMOGRAPHY (CT), AND CT DEVICE 



DECLARATION UNDER 37 C.F.R. § 1.131 



COMMISSIONER OF PATENTS 
P.O. BOX 1450 

ALEXANDRIA, VA 22313-1450 
Sir 



We, Herbert Bruder, Thomas Flohr, Annabella Rauscher-Scheibe, Karl 
Schwarz and Karl Stierstorfer, declare as toventors of the subject matter 

claimed in the above-identified application that: 



1. This Declaration is submitted as evidence that the subject 
matter disclosed and claimed in U.S. Patent Application No. 10/667,475 
was invented by us prior to July 10, 2002, which is the effective filing date of 
U.S. Patent No. 7,245,755 Bl ("Pari") entitled "Algorithm for Image 
Reconstruction and Image Noise Analysis in Computed Tomography." 

2. We conceived, and are the named inventors of, the subject 
matter disclosed and originally claimed in U.S. Patent Application No. 
10/667,475. 



Application No.: 10/667,475 
Attorney Docket No.: 32860-000622/ US 

3. U.S. Patent Application No. 10/667,475 is currently assigned to 
Siemens Aktiengesellschaft (hereinafter interchangeably referred to as "the 
company") as evidenced by the assignment recorded on Reel 014973, Frame 
0051 in the U.S. Patent and Trademark Office. 

4. The company is a company doing business in the Germany, a 
World Trade Organization (WTO) member country. 

5. Prior to July 10, 2002, we conceived a method of creating 
images in computer tomography (CT) and a CT device in Germany, which at 
the time of conception was a WTO member country. (A copy of an internal 
Invention Disclosure Form and English translation thereof are attached as 
Exhibit A). 

6. The Invention Disclosure Form was submitted on April 30, 
2002, handed over to a supervisor. Dr. Stephan Schaller. on May 7, 2002, 
and was received in the Patent Department for the company on May 23, 
2002. 

7. The subject matter disclosed and claimed in U.S. Patent 
Application No, 10/667,475 was disclosed in the Invention Disclosure Form 
entitled •Gewichtungsverfahren fur 3D-Ruckprojeckaonsverfabxen fur die 
Mehrzeller-Spiral-CT" shown in Exhibit A. 



Application No.: 10/667,475 
Attorney Docket No.: 32860r000622/US 

8. Prior to July 10, 2002, outside counsel Pr. Arniin Binder) was 
requested to prepare a patent application based on our invention (A copy of 
the order letter dated July 3, 2002 and an English translation thereof are 
attached as Exhibit B). 

9. On August 28, 2002. Dr. Binder forwarded the draft patent 
application to one of us. Dr. Karl Stierstorfer, for review (A copy of a letter 
from Dr. Binder to Dr. Stierstorfer and an English translation thereof are 
attached as Exhibit C). 

10. After reviewing the draft application, comments regarding the 
draft application were sent to Dr. Binder (Copies of emails dated August 29, 
2002 and September 4. 2002 and English translations thereof are attached 
as Exhibits D and E, respectively). 

11. On September 23. 2002, German Patent Application 10244181.2 
(the "'181 Application), which conformed substantially to the draft forwarded 
from Dr. Binder to Dr. Stierstorfer, but with the comments shown in 
Exhibits D and E, was filed with the German Patent Office (A copy of the 
'932 Application and an English translation thereof are attached as Exhibit 

n 

12. On September 23. 2003, the above-identified application was 
filed as U.S. Patent Application No. 10/667.475. and claimed priority under 
35 U.S.C § 119 to the '181 Application, The U.S. application is an English 



Application No.: 10/667.47S 
Attorney Docket No.: 32860-000622/US 

translation of the filed German Patent Application, slightly modified to place 

it in better form for U.S. practice. 

13. We hereby declare that all statements made hereto of our own 
knowledge are true and that all statements made on information and belief 
are believed to be true; and further that these statements were made with 
the knowledge that willful false statements and the like so made are 
punishable by fine and imprisonment, or both, under Title 18, United States 
Code, § 1001 and that such willful false statements may jeopardize the 
validity of the application or any patent issued thereon. 
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Dated 
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EXHIBIT A 



Vertrauljch! , 

An 

Siemens AG 

bzw. Beteiligungsgesellschaft 



ERFINDUNGSMELDUNG 

Bitte verschlossen weitersenden! 

Bereits vorab an ZT PA Qbermittelt per FAX □ 
Wennja - bitte unbedingt ankreuzen! 



Aktenzeichen der PA 



2: a 



Ich/Wir (Vor- und Nachname der/des Erfinder(s] - weitere Angaben und Unterschrift(en] letzte Seite) 

Dr. Karl Stierstorfer, Dr. Annabella Rauscher, Dr. Herbert Bruder, Dr. 
4Karl Schwarz, Dr. Thomas Flohr 



Anzahl der 
Ertinder 



Datum der Ausfertigung: 



30.4.2002 



melde[n] hiermit die auf den folgenden Seiten vollstandig beschriebene Erfindung mit der Bezeichnung: 
Gewiohtungsverfahren fur 3D-Ruckprojektionsverfahren fur die Mehrzeiler-Spiral-CT 



c/ > 
cr I 



An Vorgesetzten der/des Erfinder[s] 

Herrn/Frau Dr. Schaller 



CTC 



mit der Bitte, die nachstehenden Fragen zu beantworteri: 
a) Wann ging die Erfindungsmeldung bei Ihnen ein? 



(Dienststelle) 



b) Geht die Erfindung auf offentlich geforderte Arbeiten zuruck? 
[Xlnein | | ja, Vorhaben: 



c) Gibt es ein zugeh6riges internes FuE-Projekt? 
nein Q ja, Projekt: 



Eingang am: 



Dr. Stefan Schauta 

Ab Eingang ISuft gesetzliche Frist! 



Nur bei ZT-Erfindungen auszufullen: 

Projekt-Nr. Titel: 



Kerntechnologie: 



Entwicklungs- 
projekt 

Forschungs- 
projekt 



im Interesse von Bereich: 



Ansprechpartner 



d) Anmeldung wird empfohlen □ nein 



Kosten tragt (Organisationseinheit): 



Drlnglichkeitsvermerk 



I I Die Erfindung betrifft nicht unser Interessengebiet. Es sind noch folgende 
Dienststellen zu befragen: 



(Datum) 



(Unterschrift des Vorgesetete n) 



Bitte wegen gesetzlicher Frist sofort weiterleiten 



An 



Eingang am: 



ZT PA (Patentabteilung) 



Standort:. 



(z.B.: Mch/M, Erl/S, Bln/N, Khe/R) 

zur weiteren Veranlassung. 



CT IPS AE 



Eingang 23. Mai 2002 



GR 



1. Technisches Problem, das durch die Erfindung gelost wird 

Bei der Spiral-CT mit groBen Kegelwinkeln kflnnen Datenredundanzen durch ein Gewichtungsverfahren 
berUcksichtigt werden (vgl. Erfmdungsmeldung 2001P22782DE, Anmeldung eingereicht). Dabei treten 
jedoch charakteristische Artefakte auf. 

2. Bisherige Losung 

Fur die Behandlung von Datenredundanzen bei der Kegelstrahl-Spiral-CT sind in der Literatur keine 
Verfahren verfubbar. Bei moderatem Kegelwinkel kann auf Verfahren wie AMPR (Schaller et al., Novel 
approximate approach for high-quality image reconstruction in helical cone beam CT at arbitrary pitch 
Proc. SPIE Medical Imaging Conf. 4322) Oder SMPR (Proc. of the sixth international meeting on fully 
three-dimensional image reconstruction in radiology and nuclear medicine - 3D-2001, Pacific Grove, 2001, 
to be published) zuruckgegriffen werden. 

3. Erfindung 

Das Verfahren beruht auf einer Erweiterung des Verfahrens in Erfmdungsmeldung 2001P22782DE, wobei 
zusatzlich eine Gewichtung der ruckprojizierten Daten eingefuhrt wird: 

DieDaten P(<X,f$,q) werden nach dem Parallelrebinning nach P($,fi,q) oder P{0^p 7 q) in 
Richtung der Detektorzeilen, also in /? - bzw. p -Richtung gefiltert. Dabei ist 

(X der Fokuswinkel, 

f} der Facherwinkel, 

q die z-Koordinate (Zeilennummer), 

0=a + fi der Parallelfacherwinkel 

p= R F Sin(/F) die Parallelkoordinate (Abstand des Strahls von der 

Drehachse) und 
R F der Fokusbahnradius. 

Das entspricht einer Filterung des durch das Rebinning entstandenen virtuellen Paralleldetektors in 
Richtung der Spiraltangente. 

Fur jedes Voxel (x,y^)suche man fur jedes 0e[O,7r[ dieStrahlen P{6 + k7l>p^q) bzw. 

P{9 + k7t,P,q) , (ganze Zahlen k) deren Projektion entlang der Drehachse durch (x,y) geht und bilde die 
Summe 

wobei h eine geeignete Gewichtungsfiinktion, die die Schichtdicke bestimmt ist und d der Abstand (oder 
eine ahnliche Funktion) des Strahls vom Voxel (x,y,z) ist. Die Funktion ist dabei eine Funktion des 
Projektionswinkels. Sie sollte typischerweise am Anfang des 0-Intervalls, in dem das Voxel im 
Strahlengang liegt, klein sein und danach mit einer glatten Funktion (etwa cos 2 ) auf 1 ansteigen. Analog 
sollte es zum Ende des Intervalls, in dem das Voxel im Strahlengang liegt, auf Null fallen. 

h = £ £ w ( 0 + ( 0+k *. | ~ J, q )) 

ist die Gewichtssumme. Durch die Normierung mit der Gewichtssumme kGnnen beliebige Redundanzen, 
wie sie bei kleinen Pitchwerten auftreten, berUcksichtigt werden. 

SchlieBlich werden die Projektionswerte fiir die verschiedenen 0 Uber das Interval! 6 G [0,7l[ im 
betreffenden Voxel summiert. 




Altcrnativ kann das Verfahrcn auch wie folgt durchgefiihrt werden: 

P w (*) = jfLYF{iK* w V + tor, J jj J, q))P{6 + tor, || J, 0 , 

wobei FTeine Funktion der Zeilennummer ist. Sie sollte am oberen und unteren Rand (d.h. fur das maximale 
und minimale q) klein sein und dazwischen mit einer glatten Funktion (etwa cos 2 ) auf 1 ansteigen. Die 
Gewichtssumme ist dann 

H = E 2 W(qW*„ V + **, k }, q)) . 



Confidential! 

To Siemens AG 
or holding company 


Inventor Declaration 

please forward ahead of time 
sent ahead of time to the ZT PA 
via FAX [ ] 
if so, please mark! 


File No: 

2002E08528 DE 


I/we (first/last name of inventors) 
(further information and signatures) last 
page) 

Dr. Karl Stierstorfer, Dr. 
Annabella Rauscher, Dr. Herbert 
Bruder, Dr. Karl Schwarz, Dr. 
Thomas Flohr 


Number of inventors: 
5 


Date completed" 
04.30.2002 


Herewith report the invention described fully in the following pages: 




Weighting Method for 3-D Back Projection Method for Multi-Row Spiral CT 


To the section head(s) for the inventor(s): 

Mr/Ms. Dr. Schaller CTC 

(office) 

with the request to respond to the following questions 

a) when were you first notified of the invention? 

b) is the invention based on public funding? 
\x\ no fl yes, project: 

c) is there an associated internal FuE project? 
fxl no fl yes, project 


Received on: 
05.07.02 

Dr. Stefan Schaller 

legal period starts when received 


To be filled out only for ZT inventions: 

Project No Title: nuclear technology: 


fl development project region of interest: contact person 


[] research project 






d) application filing recommended []no [x] yes urgency remark 
Cost carrier: (organizational unit): 


[] The invention does not fall into our area of interest. The following offices must be questioned. 


05.07.02 


/s/ Schaller 




(date) 


(signature of section head) 




Please forward immediately because of the legal deadline 




To: ZT PA (patent department) 




Received on: 


Location: 

(e.g. Mch/M Erl/S Bln/N, Khe/R) 


CT IPS AE 
(May 23,2002) 

GR 


to be processed further 







1. 



Technical problem solved with the invention 



With spiral CTs having large conical angles, the data redundancies can be taken into consideration with a weighting 
method (see inventor declaration 2001P22782DE; application filed). The following artifacts occur in the process. 



2. Prior art solution 



No methods are available in literature for solving the data redundancies occurring with the conical-beam spiral CT. 
With a moderate conical angle, we have recourse to methods such as AMPR (Schaller et al., Novel approximate 
approach for high-quality image reconstruction in helical cone beam CT at arbitrary pitch Proc. SPIE Medical Imaging 
Conf. 4322) or SMPR (Proc. of the sixth international meeting on fully three-dimensional image reconstruction in 
radiology and nuclear medicine - 3D-2001 , Pacific Grove, 2001, to be published). 



3. Invention 

The method is based on an expansion of the method disclosed in the inventor declaration 2001P22782DE, wherein 
additionally a weighting of the back-projected data is introduced: 

The data P(a, P, q) are filtered after the parallel binning based on P (0, P, q) or (9, p, 1) in the direction of the detector 
rows, meaning in P and/or p direction, wherein: 

a the focusing angle 

P the fan angle 

q the z-coordinate (row number) 

0 = a + p the parallel fan angle 

p = R F sine(P) - the parallel coordinate (distance of beam to the axis of rotation) and 
R F the focal path radius 



This corresponds to a filtering of the virtual parallel detector that developed in the direction of the spiral 
antenna. 

For each voxel (x, y, z) we have to find 6e[0,7i[ the beams P(0+k7i, p , q) and/or 

P(9+k7i, P , q) (whole numbers k) for which the projection along the rotational axis runs through (x, y) 
and form the sum: 



P w (9) = ^Y.W{e + k7u)h{d xyi {0 + k7t 

H k q 



p 



q))P{G + k7r 



VP 



wherein h is a suitable weighting function which determines the layer thickness and d is the distance (or a 
similar function) of the beam from the voxel (x, y, z). The function W is a function of the projection angle. 
Typically it should be small at the start of the 0 interval where the voxel is in the beam path and later on 
should increase with a smooth function (approx. cosine 2 ) to 1 . Analog thereto, it should fall to zero at the 
end of the interval in which the voxel is in the beam path. 



H= EE W{6 + k7i)h(d xyz (0 + kn 



P 



k q 



VP 



q)) 



is the weighting sum. By standardizing the weighting sum, any optional redundancies such as occur with 
small pitch values, can be taken into consideration. 



Finally, the projection values for the different 9 can be summed up via the interval 9e[0,7t[ in the respective 
voxel. 



EXHIBIT B 



SIEMENS 



Binder Armin Dr. 
Neue Bahnhofstr. 16 

89335 Ichenhausen 



GERMANY 



Name Dr. Norbert Moritz 

Abteilung CT IP MED 

Telefon +49 9131 7 32556 

Fax +49 9131 7 32226 

E-Mail norbert.moritz@siemens.com 

Ihr Schreiben 

Unser Zeichen 200210627 
MTZ / MTZ 

Datum 03.07.2002 



AUFTR AG 

Erfindung 2002E08528 mit der Bezeichnung: Gewichtungsverfahren fur 3D- 
Ruckprojektionsverfahren fur die Mehrzeiler-Spiral-CT 



Sehr geehrter Herr Patentanwait, 

hiermit erhalten Sie Unterlagen zur oben richer bezeichneten Erfindung. Wir beauftragen Sie, 
eine Patentanmeldung bis zum 

30.09.2002 



auszuarbeiten. 



Bitte prufen Sie, ob diese Erfindungsmeldung zusammen mit einem Oder mehreren unserer 
separaten Auftrage mit Zeichen 2002 10630, 2002 10631 und 2002 10626 einer gemeinsamen 
Anmeldefassung zugefuhrt werden kann. 

Bitte senden Sie das beigefugte Doppel mit Ihrer Empfangsbestatigung umgehend an 
uns zuriick. 

Mit/freundlichem GruB Anlagen: 

• Empfangsbestatigung 

• Erfindungsmeldung in Kopie 

• Kopie der Mitteilung an Erfinder 

• Anmeldefassung 2001 P22782DE in Kopie 




Corporate Technology 

Corporate Intellectual Property and Functions 
Leitung: 

Dr. Horst Fischer 



Briefadresse: 
Siemens AG 
CT IP MED 
Postfach 32 20 
91050 Eriangen 



Hausadresse: 
Paul-Gossen-Str. 100 
91052 Eriangen 



Siemens Aki/engeselischatt • Voreltzendor dos Aulsfchtsrats; Kail -Hermann Baumann • Vorstand: Helnncn v. Pierer, Vorsitzender • Mirglleden Voflcer Jung, 
Edward G. Knjbasik. Rudi Umprecm. Heinz- Joacnim Netrburger. Peter Pribifla. Jurgen Radomski. Erich R, Reinhardl. Uriel J. Sha/el, Ctaus Weyrich. 
Klaus Wucherer 

Sl2 der Gesellschaft Berlin und MGnchen • Regtsiergerichl: Berfin-Chartottenburg, HRB 12300: MOnchon. HRB 6684 



Dr. Arming Binder 
Neue Bahnhofstr. 16 
89335 Ichenhausen 



GERMANY 
ORDER 

Invention 2002E08528 with the title: "Weighting Method for 3D-Back Projection 
Method for the Multi-Line Spiral CT 

Dear Attorney: . 

Enclosed are the documents for the above-referenced invention. We herewith 
authorize you to prepare a patent application by 

09.30.2002. 

Please check whether this declared invention can be combined into a single 
application with one or several of our separate orders with the File Numbers 2002 10630; 
2002 10631 and 2002 10626. 

Please return the enclosed second copy immediately along with your 
acknowledgment of receipt. 

Sincerely, Enclosures 

acknowledgment 

copy of invention disclosure 

copy of notice to inventor 

copy of application document for 

2001P22782DE 



EXHIBIT C 



Telefax 



Patentanwalt 

European Patent & Trademark Attorney 

DR. PHIL. NAT. ARMIN BINDER 

DiPL PHYS. 

Neue Bahnhofstr. 16 
D - 89335 Ichenhausen 
Tel: ++ 49(0)8223/6213 
Fax: ++ 49 (0)8223/6316 
e-mail:pa.binder@t-online.de 

An: Herrn Datum: 28. August 2002 

Dr. Karl Stierstorfer 

Telefax: 09191-18-9996 Seiten: 35 

Betreff: Patentanmeldung 

Ihr Zeichen: 200210627 / 2002 E 08528 DE 
Unsere Akte: 2002.0528 

Sehr geehrter Herr Dr. Stierstorfer. 

nachfolgend erhalten Sie den Entwurf zur oben genannten Patentanmeldung mit der Bitte 
Durchsicht. Bitte beachten Sie hierbei auch das "Merkblatt fur Erfinder'. 

Sollte Ihre Durchsicht nicht innerhalb einer Woche durchfuhrbar sein, so bitte ich urn 
umgehende Nachricht. 

Falls Sie mit dem Entwurf einverstanden sind. bitte ich um Rucksendung der 
Einvcrstandniserklarung, gegebenenfalis einschlieBlich der gewunschten Anderungen. 

Fur Ruckfragen stehe ich jederzeit gerne zur Verfugung. 



Mit freundlichen GruBen 




Anlagen 

Merkblatt fur Erfinder 
Anmeldeentwiirf 



jpucorc i merit, toic 



2EIT 

RUFNUMMER 
NPME 



28.08.02 10:08 

+49-8223-6316 

PA DR. BINDER ICHENH. 



NR. DPTE1 DATUM ZEIT 



BAUER 



ANZ PN 



GBT.-NR. MODE 



STRTUS 



813 F.7 28.08 09:54 13/06 35 +49 9191189996 



EC 08 M OK 



Telefax 



Patentanwalt 

European Patent & Trademark attorney 

DR. PHIL. NAT. ARMIN BINDER 

Dipl. Phys. 

Neue Bahnhofsir. 1 6 
D - 89335 Ichenhausen 
Tel; ++49 (0) 8223/6213 
Fax: ++49 (0) 8223/6316 
e-mail:pa.binder@t-online.de 

An: Herra Datum: 28. August 2002 

Dr. Karl Stierstorfer 

Telefax: 09191-18-9996 Seiten; 35 

Betreff: Patentanmeidung 

Ihr Zeichen: 200210627 / 2002 E 08528 DE 
Unsere Akte: 2002.0528 

Sehr geehrter Hen Dr. Stierstorfer, 



nachfolgend erhalten Sie den Entwurf zur oben genannten Patentanmeidung mit der Bitte um 
Durchsicht. Bitte beachten Sie hierbei auch das "Merkblatt fur Erfinder". 

Sollte Ihre Durchsicht rticht innerhalb einer Woche durchfiihrbar sein, so bitte ich um 
umgehende Nachricht. 

Falls Sie mit dem Entwurf einverstanden sind. bine ich um Rtlcksendung der 
Einverstandniserklarung, gegebenenfalls einschlielilich der gewiinschten Anderungen. 

Fur Ruckfiragen stehe ich jeder2eit gerne zur Verfilgung. 



Mit freundlichen GriiBen 




Anlagen 

Merkblatt fur Erfinder 



Fax from Dr. Phil. Nat. Armin Binder 



To: Dr. Karl Stierstorfer Date: August 28, 2002 

Fax:09191-18-9996 Pages: 35 

Reference: Patent Application 

Your File: 200210627/2002E 08528 DE 
Our File: 2002.0528 

Dear Dr. Stierstorfer: 

Attached is the draft for the aforementioned patent application with the request to 
review this draft. Please also note the "Instruction for Inventors." 

Please let me know immediately if you cannot review this text within one week. 

If you agree with the draft, please return the Consent Declaration, if applicable 
along with the desired changes. 

I will be happy to respond to any further questions. 
Sincerely, 

/s/ Dr. Armin Binder 
Enclosures: 

Instructions for inventors 
Application draft 



EXHIBIT D 



mailbox:///C|/Dolcumentc%20und%20Einstellungen/binder/Anwc... 

• • 

Subject: Anmeldung 2002.0528 

From: Sticrstorfcr Karl <kari.sticrstorfer@siemens.com> 

Date: Thu, 29 Aug 2002 12:32:46 +0200 

To: m pa.binder@t-online.dc m <pa.binder@t-on!ine.de> 

X-UIDL: I21bafe46e39746e 
X-Mozilla-Status: 0001 
X-MoziIIa«.Status2: 00000000 
Return-Path: <karl.stierstorfer@siemens.com> 

Received: from atlantis.erlm.siemens.de ([212. 1 14.202.14J) by mailin04.sul.t-online.de with csmtp id 17kMbY-01EHoWC; 
Thu, 29 Aug 2002 12:32:48 +0200 

Received: from tcenter.erlm.siemcns.de (tcenter.erlm.siemens.de [146.254.160,14]) by atlantis.erlm.siemens.de 
(8.9.3/8.9.3) with ESMTP id MAA18975 for <pa.binder@t-online.de>; Thu, 29 Aug 2002 12:32:48 +0200 (MET DST) 
Received: from erlh910a.mcd.siemens.de (ERLH910A.med.siemens.de [146.254.61.193]) by tcentcr.erlm.siemens.de 
(8.9.0/8.9.0) with ESMTP id MA A 10898 for <pa.binder@t-online.de>; Thu, 29 Aug 2002 12:32:47 +0200 (MET DST) 
Received: by erlh910a.med.siemens.de with Internet Mail Service (5.5.2656.59) id <R6KCLP36>; Thu, 29 Aug 2002 
12:32:50 +0200 

Message-ID: <CB0365D517B7D61 lB5E100508B9498B60F258E@erih904a.mcd.sicmcns.de> 
MIME-Version: 1.0 

X-Mailcr: Internet Mail Service (5.5.2656.59) 
Content-Type: text/plain; charser= , 'iso-8859-P 
Content-Transfer-Encoding: quotcd-printable 



Hallo Herr Binder, 

vielen Dank fur Ihre Anmeldung, die ich gerade durchsehe. Bin paar 
Anmerkungen babe ich da z u s 

Allgemein: sollte da nicht auch die Zeitgewichtung (Cardio) mix rein? Oder 
kommt das in eine extra Anmeldung? 

S. I, Z, 17: "bekannt" zuviel 

S. 2, 2. 21 ft; da sind m.E. zwei Dinge vermischt: zum einen dOrfen 
natOrlich fur die Rekonstruktion eines voxels nur Strahlen berucksichtigt 
werden, die sehr nahe bei dem Voxel liegen (dies wird die Gewichtung mic 
h(dx,y,2) erreicht) - das muS aber jede 3D-Ruckprojektion tun. Wichtig ist 
hier, daS die Gewichtungsfunkt ion W eingefuhrt wird, die Strahlen abhangig 
von ihrer Position ira Facher {also unabhangig von der Voxelposit ion) 
gewichtet, und zwar so, daS Strahlen, die in Achsenrichtung (!) am Rand des 
Fachers liegen. geringer gewichtet werden. 
S. 3, Z. 20ff: dito. 

S. 3, Z. 25: "demnach" trifft hier m.E. nicht, da das vorher noch gar nicht 
erwahnt wurde. 

S, 4, Z. 27: w sollte expiizit als W{q) geschrieben werden, um deutlich zu 
machen, daS das wirklich vom Zeilenindex, also von der z- Posit ion des 
Strahls im Facher abhangt . 

S. 5, Z. 29ff: die Normierung ist m.E. unbedingt not wend ig - . 

S. 6, Z Bit: Wenn man W als W(theta+kPi) schreibt, muS die 
Gewichtungsfunkt ion so beschaffen sein, daft die Gewichte am Anfang und am 
Ende des fur das Voxel relevanten theta+kPi-Bereichs reduziert werden. (Das 
ist nicht einfach - ein Voxel kann auch mehrmals in den St rahlenkegel ein- *~\ m 

und austreten - deshalb wurde ich die Beschreibung mit W(q) vorziehen. vgl . *\ 
Seiten 15-17) » 

S. 15, 2. 30ff: siehe Bemerkung zu S. 6 

S. 19, Z. 4ff: siehe Bemerkung zu Seite 2. Was ist der Unterschied von Fig 4 
und S? 

S. 20, Z. 4ff: w hangt nicht vom Abstand des Strahls zum Voxel ab, sonde rn 
nur von der (z-l Position des Strahls im Facher. 

S. 22, Z. 25ff: dito 

S. 24, Z. 7ff: "peripher" ware mir zu ungenau. Vielleicht so: * W eine 
Wichtungsfunktion darstellt, die Strahlen mit betragsmSSig groSem 
Kegelwinkel schwacher gewichtet als Strahlen mit betragsma&ig kleinem 
Kegelwinkel ** . Dann mufite naturlich irgendwo der Kegelwinkel erkiart werden, 
etwa als Winkel des Strahls zu einer axialen Ebene. 
S. 27, z. 26ff: wie Bemerkung zu S. 2 



viele GruEe, 

Karl Stieretorfer 
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E-mail dated August 29, 2002 from Karl Stierstorfer to Armin Binder 
Hello Mr. Binder: 

Thank you for the application documents which I am reviewing at this time. I have a few 
comments relating to the application: 

General: Should the time weighting (cardio) not also be included? Or will that be in a 
different application? 

P. 1 , L. 1 7 "known" is too much 

P.2 5 L. 2 Iff I believe two things are mixed up here. For one thing, for the 

reconstruction of a voxel only rays should be taken into consideration which are 
very close to the voxel (achieved for the weighting with h(dx,y,z) - however that 
must be done for any 3-D back projection. It is important here that the weighting 
function W is introduced, that the rays are weighted in dependence on their 
position in the fan (meaning independent of the voxel position). That is to say, 
rays which are located at the edge of the fan as seen in axial direction (!) are 
weighted lower. 

P. 3, L. 20ff ditto 

P. 3, L. 25: I believe "accordingly" is not appropriate here since it was not previously 
mentioned. 

P.4, L.27: W should explicitly be written as W(q) to make it clear that it actually 
depends on the line index, meaning the z position of the ray in the fan. 

P.5,L.29ff: I believe the standardization is absolutely necessary. 

P.6, L. 8ff: If we write W as W(theta+kPi), the weighting must be such that the 

weights are reduced at the start and the end of the theta+kPi region that is relevant 
for the voxel. (That is not easy - a voxel can also enter or leave a beam several 
times - 1 would therefore prefer the description using W(q); see also pages 1 5- 
17). 

P. 1 5,L.30ff: See comment to page 6 

P. 19, L.4ff: See comment to page 2. What is the difference between Figures 4 and 5? 
P.20, L.4ff: W does not depend on the distance between the ray and the voxel, but only 

on (z-)position of the ray in the fan. 
P.22, L.25ff: ditto 

P. 24, L.7ff: "Peripheral" would not be accurate enough for me. Perhaps you could 

use: " A weighting function is shown with W which gives less weight to rays with 
large conical angle." Of course, the conical angle would then have to be 
explained somewhere, e.g. as angle between the ray and an axial plane. 

P. 27, L. 26ff: See comment to P. 2 



Sincerely, 

Karl Stierstorfer 
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Subject: Rc: Ihrc Akie 2002.0528 

From: "Dr. Armin Binder" <pa.binder@t-onlinc.de> 

Date: Wed. 04 Sep 2002 16:3 i :29 +0200 

To: Stierstorfer Karl <karl. stieretorfer@siemens.com> 

X-Mozilla-Status: 0001 

X-Mo/jlla-Status2: 00000000 

Message-ID: <3D761 94 1.2070504@t-oni inc.de> 

User-Agent; Mozilla/5.0 (Windows; U; Windows NT 5.0; de-DE; rv:0.9.4. 1 ) Gecko/20020508 Netscape6/6.2.3 
X-Accept-Language: de-DE 
MIME-Vcrsion: 1.0 

References: <CB0365D5 1 7B7D61 1 B5E100508B9498B60F25B4@cr!h904a.mcd,siemens,de> 
Content-Type: text/plain; charsct=lSO-8859-l; formar=flowcd 
Conlent-Transfer-Encoding: Sbit 

Hallo Herr Stierstorfer, 

ich habe die Anmeldung nochmals wegen Ihrer Anmerkung zu W{..) durchgesehen . Ich lag of fensicht lich 
falsch. da dieser Term erst hier neu eingefOhrt ist. Ich habe deshalb - mit Ausnahme der beiden 
expliziten Beispiele in den Unteranspruchen 7 und 8 - alles in W(q) geSndert. 

Die Anderungen in den Formelrs (hHd...) in h(..)} habe ich ebenfalls durchgefuhrt . 

mfg 

Binder 

Stierstorfer Karl wrote: 
Hallo Herr Binder, 

vielen Dank fur Ihre neue Annie 1 dung, die ich jetzt leserlich erhalten habe. 

Die wichtigsten Punkte sind jetzt korrekt, nur ein paar Dinge fallen mir 
noch auf: 

S. 6, Formeln: hier sollte man niche h* <dx , y , z i . . . ) } schreiben, da man sonot 
denken konnte, daS das ein Produkt aus h und d ist, sonde rn h (dx, y, z { . . . ) ) 
S . 15, Forme 1 : dito 
S. 24. Forme! : di to 
S. 25, Formeln; dito 

S. 16, Z. 3: wenn Sie w(theta>kpi) schreiben, ist das eine 
Wtchcungsfunktion, die vom Projektionswinkel theta+kpi abhangt . DaS das, wie 
Ihr Text besagt, Strahlen bezQglich ihrer z-Anordnung gewichtet, ist nicht 
selbstverstSndlich und geht aus dem Text nirgendwo hervor . Im allgemeinen 
vurde ich die Formulierung Ober W(q) vor Ziehen, da sie intuitiver ist. 

Die Einverstandniserkl5r-ung schicke ich Ihnen unterschrieben per Fax zurOck 

Viele GruSe, 

Karl stierstorfer 
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E-mail of September 4, 2002 from Armin Binder to Karl Stierstorfer 



Dear Mr. Stierstorfer: 

I have again reviewed the application because of your comment regarding W (..). I was 
obviously wrong since this term has only here been introduced new. I therefore changed 
everything to W (q), with the exception of the two explicit examples in the dependent 
claims 7 and 8. 

I also made the changes in the formulas from (h* (d. . .) to h(..)). 
/s/ Binder 



From Karl Stierstorfer: 
Dear Mr. Binder: 

thank you for the application text which I have now received in a readable form. 

The most important points are now correct, with the exception of a few things that I 
noticed: 

P. 6 formulas: should not read h* (dx, y, z(. . .)) since it would be possible to consider this 
a product of h and d; instead it should read h (dx,y,z (. . .)). 
P. 15 formula: ditto 
P. 24 formula: ditto 
P. 25 formulas: ditto 

P. 16, line 3: if you write W(theta+kpi), that refers to a weighting function which 
depends on the projection angle theta+kpi. It is not understood that this refers to the 
weighting of rays with respect to their z arrangement, as your text claims, and also does 
not follow anywhere from the text. In general, I would prefer the formula with W(q) 
because it is more intuitive. 

I will return the signed consent declaration via facsimile. 
Sincerely, 



/s/ Stierstorfer 
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Method of creating images in computer tomography, and 
CT device for carrying out the method 

The invention relates to a method of creating images in 
5 computer tomography, in which an object under 
examination is scanned with the aid of a beam 
originating from a focus and a flat detector array, the 
output data determined is filtered in a suitable way 
and backpro j ected in order to obtain at least one slice 
10 which represents the absorption values of the section 
of the object under examination. Furthermore, the 
invention relates to a CT device which is suitable for 
carrying out this method. 

15 Methods of this type have become known under the term 
"filtered backpro j ect ion" but, as 3D methods, that is 
to say in conjunction with a matrix-like detector 
array, do not furnish an image quality which is judged 
to be adequate in practice, since what are known as 

2 0 "cone beam artefacts" occur because of the conical x- 
ray beam. 

In addition, the disadvantage with this method is that 
redundant data, such as are produced during spiral 
25 scanning with low table advance as a result of multiple 
irradiation of one and the same voxel, is not used, 
which results in the radiation dose administered to the 
object under examination being used only incompletely 
for imaging. 

30 

Furthermore, there are thoughts relating to proceeding, 
in conjunction with 2D methods for image 
reconstruction, in such a way that preliminary images 
are calculated in large numbers by means of "filtered 
35 backpro j ection" from initial data which originates from 
sections of the focal path which are intrinsically 
inadequate for image reconstruction, the preliminary 
images being reformatted to form a final slice only in 
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a second step. These 2D methods are less useful for 
detector arrays with a large width, that is to say a 
great extent in the direction of the system axis, since 
an extremely large number of preliminary images then 
has to be processed, which is a problem even when high 
computing power is available. 

In all the aforementioned methods, the problem arises 
that, because of possible data redundancy which occurs 
during the scanning of the object under examination, 
image artefacts are produced which have a detrimental 
influence on the image quality. 

The invention is based on the object of specifying a 
method of the type mentioned in the introduction which 
permits the image quality to be increased. 

According to the invention, this object is achieved by 
the features of the independent patent claims . Special 
configurations and improvements emerge from the 
features of the subclaims . 

The inventors have discovered that in image processing 
it is firstly necessary to take account of the 
redundancies which occur with respect to a respectively 
considered voxel in the region under examination, in 
particular during the backpro j ection of the filtered 
data. Secondly, however, it is also necessary to take 
into account how the ray is positioned in the beam. For 
example, it has a decisively positive effect on the 
image quality if rays which are located centrally, 
based on the extent of the beam in the direction of the 
axis of rotation, are taken more into account when 
creating the image than rays which are positioned only 
marginally in relation to the beam, based on the same 
extent . 
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Accordingly, the inventors propose an improved method 
for creating images in computer tomography which has 
the following method steps: 

- to scan an object under examination with a beam 
5 originating from at least one focus and with a flat 

detector array having a large number of distributed 
detector elements for detecting the rays of the 
beam, the at least one focus is rotated relative to 
the object under examination on at least one focal 

10 path running around the object under examination 

with opposite detector array, the detector elements 
of the detector array supplying initial data which 
represents the attenuation of the rays as they pass 
through the object under examination, 

15 - the initial data is filtered, 

- the filtered initial data is then backpro j ected 
three-dimensionally in order to produce at least one 
slice of a layer of the object under examination 
which has a layer thickness, the slice representing 

20 absorption values, obtained from the initial data, 

of the voxels belonging to the layer of the object 
under examination for the radiation of the beam, and 

- during the backpro j ect ion, the rays are weighted as a 
function of their position in the beam. 

25 

In a particular embodiment, the filtering is carried 
out in the direction of the tangent to the focal path 
belonging to the respective focal position. It has been 
shown that, with this filtering direction, a 

30 particularly high image quality can be achieved. The 
selection of this filtering direction is based on the 
finding that the abovedescribed 2D method permitting a 
high image quality and based on preliminary images 
would, so to speak, change into a 3D method if the 

35 sections of the focal path on which the calculation of 
preliminary images is based were to be shortened so 
severely that they would cover only a single projection 
whose data would then be filtered in the direction of 
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the tangent to the focal path, and in that it would 
then be possible to expect that such a 3D method would 
then permit a comparatively good image quality as the 
2D method. 

5 

If the beam has an extent in the direction of rotation 
and an extent in the direction of the axis of rotation, 
then it is advantageous to weight rays which are 
arranged centrally in the beam, as based on the extent 
10 of the beam in the direction of the axis of rotation, 
to a greater extent than the rays which are arranged 
close to the edge in the beam, as based on the extent 
of the beam in the direction of the axis of rotation. 

15 The method according to the invention may then be 
carried out particularly easily if, 'before filtering, a 
conversion of the initial data obtained in fan beam 
geometry in the form of rays P(a,$,q) into parallel 
data present in parallel beam geometry in the form of 

20 rays P(d,/3,q) (azimuthal "rebinning" ) or P(6,p,q) 
(complete "rebinning", that is to say azimuthal and 
radial "rebinning" ) is performed. In this case, with 
reference to fig. 3 

25 a is the focal angle 

(3 is the fan angle 

qr is the row index of the detector system 

corresponding to the z coordinate, 

8 = a + (3 is the parallel fan angle, 

30 p = R F sin(f3) is the parallel coordinate corresponding 

to the distance of the ray from the axis 
of rotation (system axis), and 

R F is the radius of the focal path. 

35 According to a preferred embodiment of the invention, 
the backprojection of the parallel data is carried out 
in that, in the course of the backprojection for each 
voxel V(x,y,z), for each ae[o,7r] for the rays 
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P(e + kn,p,q) and P (0 + kn,p,q ) whose projection along the 
system axis goes through (x,y), the sum 



is formed, where 

x,y,z are the coordinates of the respective voxel 
V(x, y f z) , 

k is a whole number corresponding to the number 

of half revolutions of the focus included in 
the reconstruction, 

p are the parallel coordinates of those rays 

whose projections along the system axis run 
through the coordinates (x,y) of the respective 
voxel V(x, y, z) , 
P are the fan angles of those rays whose 

projections along the system axis run through 
the coordinates (x,y) of the respective voxel 
V(x,y, z) , 

h is a weighting function determining the layer 

thickness of the layer of the object under 
examination represented in the slice produced, 

d is a function which is equal to the distance of 

the respective ray from the corresponding voxel 
V(x,y) and or is dependent on the distance of 
the respective ray from the corresponding voxel 
V(x, y) , and 

W represents a weighting function which weights 

rays with a large parallel fan angle 0 less 
than rays with a small parallel fan angle ©. 

In this case, the notation |~| expresses the fact that 

the formation of the sum can optionally be carried out 
for rays obtained by means of azimuthal u rebinning" or 
complete "rebinning" , the filtering tangential to the 
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focal path in the case of azimuthal "rebinning" being 
filtering in the (3 direction and, in the case of 
complete "rebinning" , being filtering in the p 
direction. 

5 

Because of the summing both over k and over g, it is 
then ensured that all the rays running through one and 
the same voxel are taken into account and the radiation 
dose supplied to the object under examination is thus 
10 used completely. 

A particularly preferred embodiment provides that, for 
the purpose of the backpro j ection of the parallel data, 
the sum 

15 




normalized to the sum H of the weights h 

n k g 

20 

is formed. This procedure permits an image quality 
which is improved again, since possible overemphasis of 
voxels which are struck by more rays than other voxels 
2 5 is eliminated, and thus corresponding artefacts are 
avoided. The CT value of the respective voxel is 
obtained by summing over 0. 



6+kn 



\ 



According to the invention, the weighting function can 
3 0 be a function of the parallel fan angle with W(d+ht) , 
and preferably represent a smooth function which has 
the value 1 for central rays, as based on the extent of 
the fan of rays in the direction of the axis of 
rotation or z direction, and tends to 0 for ' marginal 
3 5 rays . 
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If the detector array has detector elements arranged in 
the manner of rows, then the weighting can be 
represented by a function of the row number W (g) , a 
5 smooth function also preferably being .used here which, 
for rays with a centrally placed detector row or a 
plurality of centrally placed detector rows, has a high 
value, preferably the value 1, and for rays arranged at 
the margin tends to 0; for example this can be a cos 2 
10 function. 

The method according to the invention can be used, for 
example when, according to one variant of the 
invention, the focal path is a circular path (tomogram 

15 scanning) . According to a preferred variant of the 
invention, however, the focal path is a spiral path, 
which is brought about by the focus being moved on a 
circular path about the system axis and, at the same 
time, there being a relative movement between focus and 

2 0 object under examination in the direction of the system 
axis. On the basis of such spiral scanning, even 
• relatively large volumes of the object under 
examination can be examined without problems. 

2 5 In the case of tomogram scanning, k is normally k=l or 

k=2 ; in the case of spiral scanning, k is selected 
while taking account of the relative displacement in 
the direction of the system axis carried out per full 
revolution, such that the region to be imaged of the 

3 0 object under examination is registered completely. 

According to the basic idea of the invention, the 
inventors also propose to improve a CT device for 
scanning an object under examination with a beam 
35 originating from at least one focus and with a flat 
detector array having a large number of distributed 
detector elements for detecting the rays of the beam, 
the at least one focus moving relative to the object 
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under examination on at least one focal path running 
around the object under examination with the opposite 
detector array, to the effect that at least means for 
collecting detector data, filtering and backpro j ection 
are provided, which carry out the method described 
above, these means preferably and as far as possible 
being implemented by programs or program modules. 

The invention will be explained in more detail below 
using an exemplary embodiment which is illustrated in 
the appended schematic drawings, in which: 

Fig. 1 shows a CT device having a plurality of rows of 

detector elements in a partly perspective, 

partly block-diagram illustration; 
Fig. 2 shows a longitudinal section through the device 

according to Fig. 1; 
Fig. 3 shows a diagram illustrating "rebinning"; 
Fig. 4 shows scanning with opposite focus and detector 

in the beam with a rotating gantry; 
Fig. 5 shows the course of two different weighting 

functions W; 

Fig. 6 shows scanning with opposite focus and detector 
in the beam with the gantry rotating and moving 
in the z direction. 

Figs 1 and 2 illustrate a CT device of the third 
generation suitable for carrying out the method 
according to the invention. Its measurement 
arrangement, designated overall by 1, has an x-ray 
source, designated overall by 2, with a beam aperture 
stop 3 arranged in front of the latter, close to the 
source (Fig. 2) and a detector system 5 formed as a 
two-dimensional array of a plurality of rows and 
columns of detector elements - one of these is 
designated 4 in Fig. 1 - with a beam aperture stop 6 in 
front of the latter, close to the detector (Fig. 2). In 
Fig. 1, for clarity, only eight rows of detector 
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elements 4 are shown. However, as indicated by dots in 
Fig. 2, the detector system 5 has further rows of 
detector elements 4. 

5 The x-ray source 2 with the aperture stop 3, on the one 
hand, and the detector system 5 with the aperture stop 
6, on the other hand, are fitted, in the manner which 
can be seen from Fig. 2, to a rotary frame 7, opposite 
each other in such a way that a pyramidal x-ray beam 

10 originating from the x-ray source 2 during operation of 
the CT device and masked by the adjustable aperture 
stop 3, whose marginal rays are designated 8, strikes 
the detector system 5. In this case, the aperture stop 
6 is adjusted in accordance with the cross section of 

15 the x-ray beam set by the aperture stop 3 such that 
only that region of the detector system 5 which can be 
struck directly by the x-ray beam is- exposed. In the 
operating mode illustrated in Figs 1 and 2, these are 
eight rows of detector elements 4, which are referred 

20 to as active rows in the following text. The further 
rows, indicated by dots, are covered by the aperture 
stop 6 and therefore not active. 

Each row of detector elements 4 has a number K of 
25 detector elements, fi k = /3j to (5 K being the channel index 
and each detector element being assigned a fan angle 
(3 k . The fan angle of the central detector element is 
equal to 0; the fan angles of the two outermost 
detector elements are f5 x = +/3 max and (3 K = -/3 max . 

30 

The active rows L Q of detector elements 4 are 
designated L± to L Qt q = 1 to Q being the row index 
which, in the case of the exemplary embodiment 
described, corresponds so to speak to the z coordinate. 

35 

The x-ray beam has the cone angle cp plotted in Figs 1 
and 2, which is the opening angle of the x-ray beam in 
a plane containing the system axis Z and the focus F. 



200210627 - 10 - 

The opening angle of the x-ray beam in a plane at right 
angles to the system axis Z and containing the focus F 
(fan opening angle) is 2(5 max and is plotted in Fig. 1. 

5 The rotary frame 7 can be set rotating by means of a 
drive device 22 about a system axis designated by Z. 
The system axis Z runs parallel to the z axis of a 
three-dimensional rectangular coordinate system 
illustrated in Fig. 1. 

10 

The columns of the detector system 5 likewise run in 
the direction of the z axis, while the rows, whose 
width b is measured in the direction of the z axis and 
is 1 mm, for example, run transversely with respect to 
15 the system axis Z and the z axis. 

In order to be able to move an object under 
examination, for example a patient, into the beam path 
of the x-ray beam, a mounting device 9 is provided, 

20 which can be displaced parallel to the system axis Z, 
that is to say in the direction of the z axis, to be 
specific in such a way that there is synchronization 
between the rotational movement of the rotary frame 7 
and the trans lational movement of the mounting device, 

25 with the effect that the ratio between translation and 
rotational speed is constant, this ratio being 
adjustable by a desired value for the advance v of the 
mounting device per revolution of the rotary frame 
being selected. 

30 

It is therefore possible for a volume of an object 
under examination located on the mounting device 9 to 
be examined in the course of a volume scan, it being 
possible for the volume scan to be carried out in the 
3 5 form of spiral scanning, with the effect that, with 
simultaneous rotation of the measuring unit 1 and 
translation of the mounting device 9, a large number of 
projections from different projection directions is 
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recorded by means of the measuring unit for each 



scanning, the focus F of the x-ray source moves 
relative to the mounting device 9 on a spiral path 



extend in the a direction over at least 7t^2fi jnaXf in order 
to permit the complete reconstruction of a CT image for 
each row of detector elements, but can also be longer 
as desired within the technical limits of the CT 
10 device. 

However, because of the fact that a plurality of rows 
of detector elements 4 are present, a volume of the 
object under examination can also be examined in the 

15 course of what is known as a tomogram scan, in which no 
relative movement takes place in the direction of the z 
axis between the measuring unit 1 and mounting device 9 
(v = 0). In the case of the tomogram scan, the size of 
the volume examined is therefore determined by the 

20 number of active rows of detector elements 4. During a 
tomogram scan, the focus F moves on a circular focal 
path which lies in a plane designated the mid-plane 
below. 

2 5 The tomogram scanning can be carried out in the form of 

a part revolution or in the form of a complete 
revolution, the part revolution covering a part 
revolution interval of at . least 7t+2f3 max (half a 
revolution plus the fan opening angle) , which permits 
30 complete reconstruction of a CT image, while a full 
revolution covers 27T. 

The measured data read out in parallel from the 
detector elements of each active row of the detector 

3 5 system 5 during the spiral or tomogram scanning and 

corresponding to the individual projections P(oc f p,q) in 
fan beam geometry is subjected to digital /analog 



revolution of the measuring unit 1. 



During the spiral 



5 designated S in Fig. 1. 



The spiral scanning must 
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conversion in a data conditioning unit 10, serialized 
and transmitted to an image computer 11. 



Following preconditioning of the measured data in a 
5 preconditioning unit 12 belonging to the image computer 
11, the resultant data stream passes to a slice 
reconstruction unit 13, which uses the- measured data to 
reconstruct slices of desired layers of the object 
under examination in accordance with a method according 
10 to the invention based on "filtered backpro j ection" and 
yet to be described in detail . 

The CT images are composed of pixels - (pixel = picture 
element) assembled in the manner of a matrix, the 

15 pixels being assigned to the respective image plane, 
each pixel being assigned a CT number in Hounsfield 
units (HU) and the individual pixels being displayed in 
accordance with a CT number/gray value scale in a gray 
value corresponding to their respective CT number. In 

20 this case, each pixel illustrates a voxel (voxel = 
volume element) of the layer illustrated in the CT 
image and belonging to the object under examination. 
Since, because of the multiplicity of rows of the 
detector system 5 and, if appropriate, the spiral 

25 scanning, measured data with respect to a plurality of 
layers of the object under examination is obtained, 3D 
data is available, ' which is subjected to 3D 
backpro j ection within the context of the invention. As 
the final result, 3D image data is available in the 

30 form of a three-dimensional matrix, for example with 
the axes y:,y,z, each element in the matrix 
corresponding to a voxel V(y:,y,z) and containing the 
gray value corresponding to the associated CT number. 
Those elements of the three-dimensional matrix which 

3 5 have the same x, y, or z value then in each case 
represent a planar slice of the layer of the object 
under examination corresponding to the defining x, y, 
or z value. 
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The images reconstructed by the slice reconstruction 
unit 13 are displayed on a display unit 16, for example 
a monitor, connected to the image computer 11. 

5 

The x-ray source 2, for example an x-ray tube, is 
supplied with the necessary voltages and currents, for 
example the tube voltage U, by a generator unit 17. in 
order to be able to set these to the respectively 
10 necessary values, the generator unit 17 is assigned a 
control unit 18 with keyboard 19, which permits the 
necessary adjustments . 

In addition, the other operation and control of the CT 
15 device is carried out by means of the control unit 18 
and keyboard 19, which is illustrated by the fact that 
the control unit 18 is connected to the image computer 
11 . 

20 Inter alia, the number Q of active rows of detector 
elements 4, and therefore the position of the aperture 
stops 3 and 6, can be adjusted, for which purpose the 
control unit 18 is connected to adjusting units 20 and 
21 assigned to the aperture stops 3 and 6. 

2 5 Furthermore, the rotation time r needed by the rotary 
frame 7 for a complete revolution can be adjusted, 
which is illustrated by the fact that the drive unit 22 
assigned to the rotary frame 7 is connected to the 
control unit 18 . 

30 

Although it is possible in principle to implement the 
method according to the invention in fan beam geometry 
as well, the CT device described is preferably operated 
in a mode in which the method according to the 
35 invention is implemented in parallel beam geometry. 

Accordingly, the data obtained during the scanning of 
the region of the body of the patient which is relevant 
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to the respective examination by means of spiral or 
tomogram scanning in fan beam geometry is first of all 
converted, in a manner known per se, into data in 
parallel beam geometry by means of a method generally 
5 designated "rebinning" . This conversion is based on 
resorting the data obtained in fan beam geometry in 
such a way that rays are removed from different 
projections recorded in fan beam geometry, and joined 
together to form a projection in parallel beam 
10 geometry. In parallel beam geometry, data from an 
interval with a length it is sufficient to be able to 
reconstruct a complete image. In order to obtain this 
data, data in fan beam geometry from an interval of 
length n+2f} max must nevertheless be available. 

15 

A projection in parallel beam geometry is illustrated 
in Fig. 3. Accordingly, all n parallel rays RPi to RP N 
of this projection assume the parallel fan angle 6 with 
respect to the x axis of the coordinate system which is 
20 illustrated in Fig. 3 and coincides with that according 
to Fig. 1. 

By using the parallel ray RPi illustrated by a 
continuous line in Fig. 3, the transition from fan beam 
25 to parallel been geometry will be explained below. 

The parallel ray RPi originates from a projection 
obtained in fan beam geometry for the focal position F x 
located on the focal path S . The central ray RF zi 

3 0 belonging to this projection in fan beam geometry, 
running through the axis of rotation 14 and therefore 
the z axis of the coordinate system, . is likewise 
plotted in Fig. 3. The focal position Y 1 corresponds to 
the focal angle Ofi; this is the angle formed by the x 

35 axis and the central ray RF Z1 . As compared with the 
central ray RF zi , the ray RPi has the fan angle /3. It is 
thus easy to see that, for the parallel fan angle 9 it 
is true that 6 = or + /?. 
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The beam spacing p from the axis of rotation 14 or the 
z axis, measured at right angles to the respective 
parallel ray, is given by p=R F sin((5) . 

5 

As becomes clear by using the central ray RP Z 
illustrated in a thickened line in Fig. 3 and extending 
through the axis of rotation 14 and the x axis, this 
ray is the central ray of a projection in fan beam 
10 geometry recorded in fan beam geometry for the focal 
position F z at the focal angle a z . Since it is true that 
/3=0 for the central ray of a projection recorded in fan 
beam geometry, it becomes clear that the following is 
true of the case of central rays: 

15 

Depending on whether an azimuthal or complete 
"rebinning" is carried out, the parallel projections 
are present in the form 

2 0 or in the form 
P(0,p, q) 

where 

a is the focal angle 

/3 is the fan angle 

25 q is the row index of the detection system 

corresponding to the z coordinate, 
e=a+fi is the parallel fan angle 

p=R F sin (ft) is the parallel coordinate corresponding to 
the ray spacing from the axis of rotation 
30 (system axis)-, and 

R F is the radius of the focal path. 

In an operating mode corresponding to a first 
embodiment of the method according to the invention, 
35 which can be selected by means of the keyboard 19, the 
CT device described operates on the basis of 
projections obtained by azimuthal "rebinning" . The data 
corresponding to these projections is filtered in the /3 
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direction, that is to say in each case in the direction 
of the tangent T belonging to the focal position of the 
central ray of the respective parallel projection (see 
fig. 3) , specifically by using one of the filter 
5 kernels customary in computer tomography, e.g. a Shepp- 
Logan or Ramachandran-Lakshminarayanan kernel . 

The parallel data filtered in this way is then 
backproj ected in that, in the course of the 
10 backproj ection, for each voxel V(x,y,z) for each <9e[0,/r] 
for the rays P{0 + kn^p,q) whose projection along the 
system axis goes through (x,y) , the sum 

p, (*) = E £ w(q) ■ *fa w fp + P. ?))■ p ( d + P. f ) 

k q 

15 

is formed, where 

x,y, z are the coordinates of the respective voxel 
V(x,y, z) , 

k is a whole number corresponding to the number 

20 of half revolutions of the focus included in 

the reconstruction, 
P are the fan angle as of those rays whose 

proj ections along the system axis run through 
the coordinates (x,y) of the respective voxel 

2 5 V(>:, y, z) , and 

h is a weighting function which determines the 

layer thickness of the layer of the object 
under examination illustrated in the slice 
produced, and 

30 d is a function which is equal to the distance of 

the respective ray from the corresponding voxel 
V(x,y) or is dependent on the distance of the 
respective ray from the corresponding voxel 
V(x,y), and 

3 5 W(q) represents a weighting function which weights 

rays which are positioned peripherally or at 
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least decentrally in the beam, as based , on the 
extent in the direction of the axis of 
rotation, less than rays which are arranged 
centrally. 

5 

As a result of the selected filtering direction and as 
a result of the summation both over k and over q, 
firstly cone beam artefacts are avoided, in the 
interest of a high dose utilization, all the rays 
10 running through a voxel V(x,y, z) are taken into 
account, but in addition, the positioning of ^the ray in 
the beam is also taken into account by means of 
appropriate weighting. 

15 The absorption value Vx,y,z associated with a voxel 
V(x,y,z) is obtained by means of summation over 6 over 
at least half a revolution, that is to say by forming 

6 

20 

The CT number corresponding to the respective 
absorption value is determined in a conventional way 
from the absorption value. 

25 In this case, different weighting functions h and 
different functions d can be set by means of the 
keyboard 19 . 

For example, a triangular or trapezoidal function is 
3 0 suitable as a weighting function h. 

The weighting function W is represented as a function 
of the projection angle 6. 

3 5 The function d can be the distance of the respective 
parallel ray from the voxel V(x,y,z) or, instead, for 
example the z (axial) component of this distance. 
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In a modification of the first operating mode described 
above, for the purpose of backpro j ection of the 
parallel data, the sum 

5 

6 

normalized to the sum • H of the weights h 

10 

is formed. This permits an image quality which is again 
improved, since the possible overemphasis of voxels 
which are * illuminated" , that is struck by rays, in a 

15 plurality of half revolutions is eliminated and thus 
corresponding artefacts are avoided. This redundancy 
occurs during spiral scans when the relative 
displacement taking place per full revolution of the 
measuring arrangement is so low (low pitch) that voxels 

20 are irradiated repeatedly. 

A second operating mode, which corresponds to a further 
embodiment of the method according to the invention and 
can be selected by means of the keyboard 19, differs 

25 from the first operating mode in that the CT device 
described does not operate on the basis of projections 
obtained by means of azimuthal but by means of complete 
"rebinning" . The data corresponding to these 
projections is filtered in the p direction in the case 

30 of complete "rebinning" , that is to say in each case 
likewise in the direction of the tangent T belonging to 
the focal position of the central ray of the respective 
parallel projection (see Fig. 3). 
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Accordingly, for the parallel data filtered in this way 
in the course of the backpro j ection, the sum 

p w (') = 17 Z Z • h i< y , ( 0 + **> fi> <*))• p ( & + k *> P> *) 

5 

is formed, p being the parallel coordinates of those 
rays whose projections along the system axis run 
through the coordinates (x,y) of the respective voxel 
V(x,y, z) . 

10 

In the case of the second operating mode, too, in a 
modification relating to the backpro j ection of the 
parallel data, a sum normalized to the sum H of the 
weights h 

15 

' p w (') = Z Z w <a) ' h b*#* & + k ^ p> *)) - p ( & + p> * ) 

k q 

is formed, namely the sum 

h = ZZ^te)* *fc,>, + k *> p> 9)) 

k q 

20 

In the case of the first and second operating modes 
described above, a function of the method according to 
the invention is provided in which, in conjunction with 
25 a voxel V(x, y, z) , all the rays whose projection along 
the axis of rotation of 14 or the z axis goes through 
x,y are considered. Whether and to what extent these 
rays are taken into account is determined by the 
weighting functions W and h and the function d. 

30 
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However, the CT device can also have further operating 
modes which can be selected by means of the keyboard 19 
and which correspond to those described previously, 
with the difference that, for a given focal position, 
5 the theoretical ray running through the respective 
voxel V(x,y,z) is determined and then, taking into 
account the weighting function h and the function d in 
the sum formation, only those rays which can actually 
make a contribution to the sum are included in'- the 
10 course of the backpro j ection . 

In the exemplary embodiments described, the relative 
movement between the measuring unit 1 and the mounting 
device 9 is in each case produced by the mounting 
15 device 9 being displaced. However, within the scope of 
the invention, there is also the possibility of 
arranging for the mounting position 9 to be fixed in 
position and, instead, to displace the measuring unit 
1. In addition, within the scope of the invention, 

2 0 there is the possibility of producing the necessary 

relative movement by displacing both the measuring unit 
1 and the mounting device 9 . 

In Figure 4, the problem of weighting the rays in a 2D 
25 representation is shown, simplified and illustrated 
schematically. In this example, the gantry of the CT 
with focus and detector in the 0° position is shown by 
continuous lines and in the 180° position, rotated 
about the z axis by 180°, is shown by dots. The 
30 detector has 12 detector rows Ll to L12, which are each 
struck by the associated rays SI to S12 from the beam. 
It should be pointed out that the spatial dimensions 
shown are illustrated exaggerated with respect to the 
extent of the detector, in order to be able to explain 

3 5 the invention more clearly. The gantry shown moves on 

a circular path here about the z axis or system axis. 
Movement in the direction of the z axis does not take 
place . 
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If the voxel V2 which lies centrally in the beam is 
considered, this voxel will be penetrated by the 
central rays S4 in both positions of the gantry. 
According to the invention, the measurements from these 
rays are highly weighted. At the same time, in the 0° 
position, the voxel VI which lies somewhat more 
peripherally is penetrated in the 0° position by the 
outer ray SI and by the ray S4 lying further in. In 
accordance with the weighting according to the 
invention, the ray SI will be given less weight than 
the ray S4 . If the likewise peripherally arranged voxel 
V3 is considered, then the weighting here is reversed. 
Since, here, the detector is constructed with detector 
elements which are divided up into rows and columns in 
the manner of a matrix, the weighting of the rays can 
be coupled directly to the row number q. 

In figure 5, located below, three examples of a course 
according to the invention of weightings W q are 
illustrated. Common to all is the greater weighting of 
the central rays with respect to the outer rays. 
However, the functions are different. The dotted line 
represents a course which drops off linearly on either 
side of the centre, the dashed line results in 
identical weighting over a large central region with a 
tendency to drop off rapidly at the periphery, and the 
continuous line shows a bell-like course which permits 
a gentle transition between the extreme values 0 and 1. 

Figure 6 shows CT beams from a gantry of a spiral CT 
respectively in the 0° position and adjacent 180° 
position. In accordance with the essence of the spiral 
CT, these positions are offset in relation to each 
other in addition to the rotation in the z direction. 
In this illustration, the voxel VI is penetrated by the 
ray SI both in the 0° position and in the 180° 
position. In accordance with this decentral position 
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of the ray SI and the associated detector row Ll 
located on the outside in the detector, these two 
measurements are weighted lowly. By contrast, the 
voxel V2 in the 0° position is penetrated by the ray 
S12 located decentrally, which is weighted lowly in 
this position for the voxel V2 , but in the 180° 
position is penetrated by the central ray S6, which is 
weighted highly for this voxel V2 . In accordance with 
the position of the rays Si in the 0° and 180° 
positions, these rays are both weighted lowly in 
relation to the voxel VI. Overall, as a result of this 
additional weighting W illustrated here, the result is 
a substantial improvement in the image quality. 

It should be pointed out that the conical x-ray beam 
has a rectangular cross section in the exemplary 
embodiment described here, but other cross-sectional 
geometries are also possible within the scope of the 
invention . 

In connection with the exemplary embodiments described 
above, CT devices of the third generation are used, 
that is to say the x-ray source and the detector system 
are displaced jointly about the system axis during the 
image generation. However, the invention can also be 
used in connection with CT devices of the fourth 
generation, in which only the x-ray source is displaced 
about the system axis and interacts with a stationary 
detector ring, if the detector system is a multi-row 
array of detector elements. 

The method according to the invention can also be used 
in CT devices of the fifth generation, that is to say 
CT devices in which the X-radiation originates from a 
plurality of foci of one or more x-ray sources 
displaced about the system axis, and not just from one 
focus, if the detector system has a multi-row array of 
detector elements . 
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The CT devices used in connection with the exemplary 
embodiments described above have a detector system of 
detector elements arranged in the manner of an 
5 orthogonal matrix. However, the invention can also be 
used in connection with CT devices whose detector 
system has a flat array with detector elements arranged 
in another way. 

10 The exemplary embodiments described above relate to the 
medical application of the method according to the 
invention. However, the invention can also be applied 
outside medicine, for example in checking luggage or in 
material examination . 
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Patent claims 

1 . A method of creating images in computer 
tomography, comprising the method steps: 

1.1 to scan an object under examination with a beam 
originating from at least one focus and with a 
flat detector array having a large number of 
distributed detector elements for detecting the 
rays of the beam, the at least one focus is 
rotated relative to the object under examination 
on at least one focal path running around the 
object under examination with opposite detector 
array, the detector elements of the detector array 
supplying initial data which represents the 
attenuation of the rays as they pass through the 
object under examination, 

1.2 the initial data is filtered, 

1.3 the filtered initial data is then backpro j ected 
three-dimensionally in order to produce at least 
one slice of a layer of the object under 
examination which has a layer thickness, the slice 
representing absorption values, obtained from the 
initial data, of the voxels belonging to the layer 
of the object under examination for the radiation 
of the beam, and 

1.4 during the backpro j ection, the rays are weighted 
as a function of their position in the beam. 

2. The method as claimed in claim 1, characterized in 
that the filtering is carried out in the direction 
of the tangent to the focal path belonging to the 
respective focal position. 

3. The method as claimed in either of claims 1 and 2, 
characterized in that the beam has an extent in 
the direction of rotation and an extent in the 
direction of the axis of rotation, rays which are 
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arranged centrally in the beam, as based on the 
extent of the beam in the direction of the axis of 
rotation, being weighted to a greater extent than 
the rays which are arranged close to the edge in 
5 the beam, as based on the extent of the beam in 

the direction of the axis of rotation. 

4. The method as claimed in one of claims 1 to 3 , 
characterized in that, before filtering, a 

10 conversion of the initial data obtained in fan 

beam geometry in . the form of rays P(a,j3,q) into 
parallel data present in parallel beam geometry in 
the form of rays P(d,f3,q) or P(6,p,q) is 
performed, where 

15 a is the focal angle 

/3 is the fan angle 

gr is the row index of the detector 

system corresponding to the z 
coordinate, 

20 6 ~ a + (3 is the parallel fan angle, 

p = R F sin((3) is the parallel coordinate 

corresponding to the distance of the 
ray from the axis of rotation (system 
axis ) , and 

25 R F is the radius of the focal path. 

5. The method as claimed in claim 4, characterized in 
that the backpro j ection of the parallel data is 
carried out in that, in the course of the 

30 backpro j ection for each voxel V(x,y,z), for each 0 

e[o,7i] for the rays P(9 + k7t,p,q ) and P(0 + kn,p,q) 
whose projection along the system axis goes 
through (x,y) , the sum 



35 
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is formed, where 

x,y,z are the coordinates of the respective voxel 
V(x,y,z), 

k is a whole number corresponding to the 

5 number of half revolutions of the focus 

included in the reconstruction, 
p are the parallel coordinates of those rays 

whose projections along the system axis run 
through the coordinates fx, y) of the 

10 respective voxel V(x, y, z) , 

j3 are the fan angles of those rays whose 
projections along the system axis run 
through . the coordinates (x,y) of the 
respective voxel V(x,y,z), 

15 h is a weighting function determining the 

layer thickness of the layer of the object 
under examination represented in the slice 
produced, 

d is a function which is equal to the distance 

2 0 of the respective ray from the corresponding 

voxel \f(x,y) or is dependent on the distance 
of the respective ray from the corresponding 
voxel V (x, y) , and 
W represents a weighting function which 

25 weights rays with a large parallel fan angle 

d less than rays with a small parallel fan 
angle 6 

6. The method as claimed in claim 5, characterized in 

3 0 that, during the backpro j ect ion of the parallel 

data, the sum 




35 



normalized to the sum H of the weights • h 
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is formed. 

5 7. The method as claimed in either of claims 5 and 6, 
characterized in that the weighting function 
represents a function of the parallel fan angle 
with W(6+hz) . 

10 -8. The method as claimed in claim 7, characterized in 
that the weighting function w(6+krr) ^represents a 
smooth function which has the value 1 for 
centrally arranged rays in the beam, based on the 
parallel fan angle, and tends to 0 for rays 

15 arranged at the edge. 

9. The method as claimed in either of claims 4 and 5, 
characterized in that the detector array has 
detector elements arranged in the manner of rows, 

2 0 and the weighting function represents a function 

of the row number W(q) . 

10. The method as claimed in claim 8, characterized in 
that the weighting function W(q) represents a 

2 5 smooth function which has the value 1 for rays to 

at least one centrally located detector row and 
tends to 0 for rays to detector rows at the edge. 

11. The method as claimed in one of claims 1 to 9, 
30 characterized in that the focal path is a circular 

path. 

12. The method as claimed in one of claims 1 to 9 , 
characterized in that the focal path is a spiral 

3 5 path which is brought about by the focus being 

moved about the system axis on a circular path 
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and, at the same time, a relative movement between 
focus and object under examination in the 
direction of the system axis taking place. 

The method as claimed in one of claims 1 to 12, 
characterized in that the detector elements on the 
detector array are arranged distributed in rows 
and lines, preferably in the manner of a matrix. 

A CT device for scanning an object under 
examination, with a beam originating from at least 
one focus and with a flat detector array having a 
large number of distributed detector elements for 
detecting the rays of the beam, the at least one 
focus moving relative. to the object under 
examination on at least one focal path running 
around the object under examination with the 
opposite detector array, at least means for 
collecting detector data, filtering and 

backpro j ection being provided as claimed in one of 
claims 1 to 13 . 

The CT device as claimed in claim 14, 
characterized in that the aforementioned 
functional means are at least partly implemented 
by programs or program modules . 
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Abstract 

Method of creating images in computer tomography, and 
CT device for carrying out the method 

The invention relates to a method for image 
reconstruction for computer tomography with a flat, 
extended detector, the rays being weighted during the 
backpro jection as a function of their position in the 
beam . 

In addition, the invention relates to a CT device for 
carrying out the above-mentioned method. 



FIG. 4 
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Beschreibung 

Verfahren zur Bilders tellung in der Computer tomographie und 
CT-Gerat zur Durchfiihrung des Verfahrens 

Die Erfindung betrifft ein Verfahren zur Bilders tellung in 
der Computer tomographie, bei dem ein Untersuchungsob j ekt mit 
Hilfe eines von einem Focus ausgehenden Strahlenbundels und 
einem flachig ausgebildeten Detektorarray abgetastet wird, 
die ermittelten Ausgangsdaten in geeigneter Weise gefiltert 
und ruckproj iziert werden, um mindestens ein Schnittbild zu 
erhalten, das die Absorptionswerte des Schnittes des 
Untersuchungsob jekts reprasentiert. Des wei.teren betrifft die 
Erfindung ein CT-Gerat, das zur Durchfiihrung dieses 
Verfahrens geeignet ist. 

Es sind Verfahren dieser Art unter dem Begriff "Filtered 
Backprojektion" bekannt geworden, die jedoch als 3D- 
Verfahren, d.h. im Zusammenhang mit einem matrixartigen 
Detektorarray, keine in der Praxis als ausreichend erachtete 
Bildqualitat erbringen, weil wegen des kegelf ormigen 
Rontgenstrahlenbiindels sogenannte "cone beam artifacts" 
auf treten. 

Aufierdem ist an diesen Verfahren nachteilig, dass redundante 
Daten, wie sie bei Spiralabtastungen mit geringem 
Tischvorschub durch mehrmaliges Bestrahlen ein und desselben 
Voxels entstehen, nicht genutzt werden, was zur Folge hat, 
dass die dem Untersuchungsob j ekt verabreichte Strahlendosis 
nur unvollstandig zur Bildgebung genutzt wird. 

Weiter gibt es Uberlegungen im Zusammenhang mit 2D-Verfahren 
zur Bildrekonstruktion so vorzugehen, dass vorlaufige Bilder 
in groSer Zahl durch "Filtered Backprojektion" aus 
Ausgangsdaten berechnet werden, die aus zur 

Bildrekonstruktion an sich nicht ausreichenden Abschnitten 
der Fokusbahn stammen, wobei die vorlaufigen Bilder erst in 
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einem zweiten Schritt zu einem endgultigen Schnittbild 
reformatiert werden. Diese 2D-Verfahren sind fur 
Detektorarrays mit groSer Breite, d.h. grofier Erstreckung in 
Richtung der Systemachse, weniger brauchbar, weil dann eine 
extrem grofie Anzahl von vorlaufigen Bildern zu verarbeiten 
ist, was selbst bei grofier zur Verfugung stehender 
Rechenleistung problematisch ist. 

Bei alien genannten Verfahren ergibt sich das Problem, dass 
aufgrund moglicher Datenredundanzen, die bei der Abtastung 
des Untersuchungsobjektes auftreten, Bildartef akte entstehen, 
welche die Bildqualitat negativ beeinf lussen . 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der 
eingangs genannten Art anzugeben, das eine Steigerung der 
Bildqualitat ermoglicht . 

Nach der Erfindung wird diese Aufgabe durch die Merkmale der 
unabhangigen Patentanspruche gelost. Besondere 
Ausgestaltungen und Verbesserungen ergeben sich aus den 
Merkmalen der Unteranspruche . 

Die Erfinder haben erkannt, dass es einerseits notwendig ist 
bei der Bildauf bereitung, insbesondere bei der Ruckpro j ektion 
der gefilterten Daten, die auftretenden Redundanzen bezuglich 
eines jeweils betrachteten Voxels im Untersuchungsbereich zu 
beriicksichtigen. Andererseits muss jedoch auch beriicksichtigt 
werden, wie der Strahl im Strahlenbundel positioniert ist. So 
wirkt es sich auf die Bildqualitat entscheidend positiv aus, 
wenn Strahlen, die zentral - bezogen auf die Ausdehnung des 
Strahlenbundels in Rotationsachsenrichtung - liegen, starker 
bei der Bilderstellung beriicksichtigt werden, als Strahlen, 
die bezogen auf die gleiche Ausdehnung nur randseitig zum 
Strahlenbundel positioniert sind. 
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Demgema.fi schlagen die Erfinder ein verbessertes Verfahren zur 
Bilderstellung in der Computertomographie vor, welches die 
folgenden Verf ahrensschritte aufweist: 

- zur Abtastung eines Untersuchungsob j ekts mit einem von 

5 mindestens einem Fokus ausgehenden Strahlenbundel und mit 

einem flachig ausgebildeten Detektorarray mit einer 
Vielzahl von verteilten Detektorelementen zum Detektieren 
der Strahlen des Strahlenbundels wird der mindestens eine 
Fokus relativ zu dem Untersuchungsob j ekt auf mindestens 

10 einer das Untersuchungsob j ekt umlaufenden Fokusbahn mit 

gegenuber 1 i egendem Detektorarray rotiert, wobei die 
Detektorelemente des Detektorarray Ausgangsdaten liefern, 
welche die Schwachung der Strahlen beim Durchgang durch 
das Untersuchungsob j ekt reprasentieren, 

15 - die Ausgangsdaten werden gefiltert, 

- die gefilterten Ausgangsdaten werden anschlieSend zur 
Erzeugung wenigstens eines Schnittbildes einer eine 
Schichtdicke aufweisenden Schicht des Untersuchungsob j ekts 
dreidimensional ruckpro j iziert , wobei das Schnittbild aus 

2 0 den Ausgangsdaten gewonnene Absorptionswerte der zu der 

Schicht des Untersuchungsob j ekts gehorigen Voxel fur die 
Strahlung des Strahlenbundels reprasentiert , und 

- bei der Ruckpro j ektion die Strahlen in Abhangigkeit von 
ihrer Position im Strahlenbundel gewichtet werden. 

25 

In einer besonderen Ausfiihrung erfolgt die Filterung in 
Richtung der zu der jeweiligen Fokusposi tion gehorigen 
Tangente an die Fokusbahn. Es hat sich gezeigt, dass mit 
dieser Filterrichtung eine besonders hohe Bildqualitat 

3 0 erzielt werden kann. Der Wahl dieser Filterrichtung liegt die 

Erkenntnis zugrunde, dass das beschriebene, eine hohe 
Bildqualitat ermoglichende, auf vorlaufigen Bildern 
basierenden 2D-Verfahren sozusagen in ein 3D-Verfahren 
iibergehen wiirde, wenn die der Berechnung vorlaufiger Bilder 
3 5 zugrundeliegenden Abschnitte der Fokusbahn so stark verkurzt 
wurden, dass sie nur noch eine einzige Pro j ektion umfassen 
wurden, deren Daten dann in Richtung der Tangente an die 
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Fokusbahn gefiltert wtirden, und dass erwartet werden konnte, 
dass ein solches 3D-Verfahren dann eine vergleichbar gute 
Bildqualitat wie das 2D-Verfahren ermoglichen wtirde . 

5 Weist das Strahlenbtindel eine Ausdehnung in Rotationsrichtung 
und eine Ausdehnung in Rotationsachsenrichtung auf, so ist es 
vorteilhaf t, Strahlen, die im Strahlenbtindel bezogen auf die 
Ausdehnung des Strahlenbundels in Rotationsachsenrichtung 
zentral angeordnet sind, starker zu wichten, als die 
10 Strahlen, die im Strahlenbtindel randnah bezogen auf die 

Ausdehnung des Strahlenbtindels in Rotationsachsenrichtung 
angeordnet sind. 

Das erf indungsgemafSe Verfahren lasst sich dann besonders 
15 leicht implementieren, wenn vor dem Filtern eine Umrechnung 

der in Facherstrahlgeometrie gewonnenen Ausgangsdaten in Form 
von Strahlen P{a,/3,q) in in Parallels trahlgeometrie 
vorliegende Paralleldaten in Form von Strahlen P(0,/3,q) 
(azimutales "rebinning") bzw. P(0,p,q) (vollstandiges 
20 "rebinning", d.h. azimutales und radiales "rebinning") 
vorgenommen wird. Dabei sind bezugnehmend auf Fig. 3 



a der Fokuswinkel 

J3 der Facherwinkel 

2 5 q der der z-Koordinate entsprechende Zeilenindex 

des Detektorsys terns # 
6 = a + f3 der Parallel facherwinkel 

p = R F sin(/?) die dem Strahlabstand von der Drehachse 

(Systemachse) entsprechende Parallelkoordinate , 

3 0 und 

R F der Radius der Fokusbahn. 



Die Rtickprojektion der Paralleldaten erfolgt nach einer 
bevorzugten Ausf tihrungsf orm der Erfindung in der Weise, dass 
3 5 im Zuge der Rtickprojektion fur jedes Voxel V(x,y,z) ftir jedes 
0e[O,/r[ ftir die Strahlen p\@ + kn,p,q) bzw. P{6 + kn,p,q) , deren 
Projektion entlang der Systemachse durch (x,y) geht, die Summe 
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k q 



d x Ao + kn\ P \q 



x>y.z 

V 




gebildet wird, wobei 
5 x,y,z die Koordinaten des jeweiligen Voxels V(x,y,z) sind, 
k eine ganze Zahl entsprechend der Anzahl der in die 

Rekonstruktion einbezogenen Halbumlaufe des Fokus, 
p die Parallelkoordinaten derjenigen Strahlen sind, deren 

Projektionen langs der Systemachse durch die 
10 Koordinaten (x,y) des jeweiligen Voxels V(x,y,z) verlaufen, 

J3 die Facherwinkel derjenigen Strahlen sind, deren 

Projektionen langs der Systemachse durch die 
Koordinaten (x,y) des jeweiligen Voxels V(x,y,z) verlaufen, 
und 

15 h eine die Schichtdicke der in dem erzeugten Schnittbild 

dargestellten Schicht des Untersuchungsobj ekts 
bestimmende Gewichtungsf unktion ist, sowie 
d eine Funktion ist, die gleich dem Abstand des 

jeweiligen Strahls von dem entsprechenden Voxel V(x,y) 

2 0 ist bzw. von dem Abstand des jeweiligen Strahls von dem 

entsprechenden Voxel V(x,y) abhangt, und 
W eine Wichtungs funktion darstellt, welche Strahlen mit 

einem groSen Parallelf acherwinkel 6 geringer wichtet 
als Strahlen mit einem kleinem Parallelf acherwinkel 0. 



Die Schreibweise \ ~l bringt dabei zum Ausdruck, dass die 



Summenbildung wahlweise fur durch azimutales "rebinning" oder 
durch vollstandiges "rebinning" gewonnene Strahlen erfolgen 
kann, wobei die Filterung tangential zur Fokusbahn im Falle 
30 des azimutalen "rebinnings" einer Filterung in /?-Richtung und 
im Falle des volls tandigen "rebinnings" einer Filterung in p- 
Richtung ist. 



35 



Infolge der Summierung sowohl iiber k als auch q ist dann 
sichergestellt , dass alle durch ein und dasselbe Voxel 
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verlaufende Strahlen beriicksichtigt werden und die dem 
Untersuchungsobj ekt zugefuhrte Strahlendosis somit 
vollstandig genutzt wird. 

Eine besonders bevorzugte Ausf uhrungsf orra sieht vor, dass zur 
Ruckprojektion der Paralleldaten die auf die Summe H der 
Gewichte h 




gebildet wird. Diese Vorgehensweise ermoglicht eine nochmals 
verbesserte Bildqualitat , da eine mogliche Uberbetonung von 
Voxel, die von mehr Strahlen als andere Voxel getroffen 
werden, beseitigt wird und somit entsprechende Artefakte 
vermieden werden. Der CT-Wert des jeweiligen Voxels wird 
durch Summation uber 0 gewonnen. 

Erf indungsgemaS kann die Wichtungsf unktion eine Funktion des 
Parallelfacherwinkels mit W(0+kn) sein, vorzugsweise eine 
glatte Funktion darstellen, die fur zentrale Strahlen - 
bezogen auf die Ausdehnung des Strahlenf achers in 
Rotationsachsenrichtung oder z-Richtung - den Wert 1 aufweist 
und fur randseitige Strahlen gegen Null geht . 

Weist das Detektorarray zeilenartig angeordnete 
Detektorelemente auf, so kann die Wichtung durch eine 
Funktion der Zeilennummer W(q) darstellt werden, wobei auch 
hier bevorzugt eine glatte Funktion verwendet wird, die fur 
Strahlen mit zentral gelegener Detektorzeile oder mehreren 
zentral gelegenen Detektorzeilen einen hohen Wert, 
vorzugsweise den Wert 1, aufweist und fur randseitig 
angeordnete Strahlen gegen Null geht, beispielsweise kann 
dies eine cos 2 -Funktion sein. 
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Das erf indungsgemafce Verfahren kann beispielsweise Anwendung 
finden, wenn gemaS einer Variante der Erfindung die Fokusbahn 
eine Kreisbahn ist ( Tomogr ammabt as tung ) . GemaS einer 
bevorzugten Variante der Erfindung ist die Fokusbahn jedoch 
5 eine Spiralbahn, die dadurch zustande kommt, dass der Fokus 
auf einer Kreisbahn urn die Systemachse bewegt wird und 
gleichzeitig eine Relativbewegung zwischen Fokus und 
Untersuchungsobjekt in Richtung der Systemachse erf olgt . Auf 
Basis einer solchen Spiralabtastung konnen problemlos auch 
10 groSere Volumina des Untersuchungsobjekts untersucht werden. 

Im Falle einer Tomogrammabt as tung gilt fur k normalerweise 
k=\ oder k=2 , im Falle einer Spiralabtastung wird k unter 
Beriicksichtigung der pro Vollumlauf erfolgenden 
15 Relatiwerschiebung in Richtung der Systemachse so gewahlt, 
dass der abzubildende Bereich des Untersuchungsobjekts 
vollstandig erfasst ist. 

Entsprechend dem Grundgedanken der Erfindung schlagen die 

2 0 Erfinder auch vor, ein CT-Gerat zur Abtastung eines 

Untersuchungsobjekts' mit einem von mindestens einem Fokus 
ausgehenden Strahlenbundel und mit einem flachig 
ausgebildeten Detektorarray mit einer Vielzahl von verteilten 
Detektorelementen zum Detektieren der Strahlen des 
25 Strahlenbiindels, wobei der mindestens eine Fokus relativ zu 
dem Untersuchungsobjekt auf mindestens einer das 
Untersuchungsobjekt umlaufenden Fokusbahn mit 
gegenuberliegendem Detektorarray bewegt, dahingehend zu 
verbessern, dass zumindest Mittel zur Sammlung von 

3 0 Detektordaten, Filterung und Ruckpro j ektion vorgesehen sind, 

die das oben beschriebene Verfahren durchfuhren, wobei diese 
Mittel bevorzugt und soweit moglich durch Programme oder 
Programm-Module verwirklicht werden. 

3 5 Die Erfindung wird nachfolgend anhand eines in den 
beigefugten schematischen Zeichnungen dargestell ten 
Ausfuhrungsbei spiels naher erlautert . Es zeigen: 
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Fig. 1 in teils perspektivischer , teils 

blockschaltbildartiger Darstellung ein mehrere 
Zeilen von Detektorelementen aufweisendes CT-Gerat; 

Fig. 2 einen Langs schnitt durch das Gerat gemafc Fig. 1 ; 

Fig. 3 ein das "rebinning" veranschaulichendes Schaubild; 

Fig. 4 Abtastung mit gegenuberliegendem Fokus und Detektor 
im Strahlenbundel bei rotierendem Gantry; 

Fig. 5 Verlauf zweier unterschiedlicher 
Wichtungsf unktionen W; 

Fig. 6 Abtastung mit gegenuberliegendem Fokus und Detektor 
im Strahlenbundel bei rotierendem und in z-Richtung 
bewegtem Gantry. 

In den Fig. 1 und 2 ist ein zur Durchfuhrung des 
erf indungsgema&en Verfahrens geeignetes CT-Gerat der 3 . 
Generation dargestellt. Dessen insgesamt mit 1 bezeichnete 
Messanordnung weist eine insgesamt mit 2 bezeichnete 
Rdntgenstrahlenquelle mit einer dieser vorgelagerten 
quellennahen Strahlenblende 3 (Fig. 2) und ein als 
f lachenhaftes Array von mehreren Zeilen und Spalten von 
Detektorelementen - eines von diesen ist in Fig. 1 mit 4 
bezeichnet - ausgebildetes Detektorsystem 5 mit einer diesem 
vorgelagerten detektornahen Strahlenblende 6 (Fig. 2) auf . In 
Fig. 1 sind der libers ichtlichkeit halber nur acht Zeilen von 
Detektorelementen 4 dargestellt. Das Detektorsystem 5 weist 
jedoch, was in der Fig. 2 punktiert angedeutet ist, weitere 
Zeilen von Detektorelementen 4 auf. 

Die Rdntgenstrahlenquelle 2 mit der Strahlenblende 3 
einerseits und das Detektorsystem 5 mit der Strahlenblende 6 
andererseits sind in aus der Fig. 2 ersichtlicher Weise an 
einem Drehrahmen 7 einander derart gegeniiberliegend 
angebracht, dass ein im Betrieb des CT-Gerats von der 
Rontgenstrahlenguelle 2 ausgehendes , durch die einstellbare 
Strahlenblende 3 eingeblendetes , pyramidenf ormiges 
Rontgenstrahlenbundel, dessen Randstrahlen mit 8 bezeichnet 
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sind, auf das Detektorsystem 5 auftrifft. Dabei ist die 
Strahlenblende 6 dem mittels der Strahlenblende 3 
eingestellten Querschnitt des Rontgenstrahlenbiindels 
entsprechend so eingestellt, dass nur derjenige Bereich des 
5 Detektorsystems 5 freigegeben ist, der von dem 

Rontgenstrahlenbiindel unmittelbar getroffen werden kann. Dies 
sind in dem in den Fig. 1 und 2 veranschaulichten 
Betriebsmodus acht Zeilen von Detektorelementen 4, die im 
Folgenden als aktive Zeilen bezeichnet werden. Die weiteren 
10 punktiert angedeuteten Zeilen sind von der Strahlenblende 6 
abgedeckt und daher nicht aktiv. 

Jede Zeile von Detektorelementen 4 weist eine Anzahl K von 
Detektorelementen auf wobei pk=P\ bis /3 K der Kanalindex ist 
15 und jedem Detektorelement ein Facherwinkel /?* zugeordnete 
ist. Der Facherwinkel des mittleren Detektorelementes ist 
gleich Null; die Facherwinkel der beiden auSersten 
Detektorelemente sind /?i=+#nax und J3 K =- #nax 

2 0 Die aktiven Zeilen L q von Detektorelementen 4 sind in Fig. 2 
mit L\ bis L Q bezeichnet, wobei q = 1 bis Q der Zeilenindex 
ist, der im Falle des beschriebenen Ausf uhrungsbeispiels 
sozusagen der z-Koordinate entspricht. 

2 5 Das Rontgenstrahlenbiindel weist den in Fig. 1 und 2 

eingetragenen Conewinkel <p auf, bei dem es sich um den 
Of fnungswinkel des Rontgenstrahlenbiindels in einer die 
Systemachse Z und den Fokus F enthaltenden Ebene handelt. Der 
Of fnungswinkel des Rontgenstrahlenbiindels in einer 

30 rechtwinklig zu der Systemachse Z liegenden und den Fokus F 
enthaltenden Ebene (Facherof fnungswinkel ) betragt 2/^ax und 
ist in Fig. 1 eingetragen. 

Der Drehrahmen 7 kann mittels einer Antriebseinrichtung 22 um 
35 eine mit Z bezeichnete Systemachse in Rotation versetzt 

werden. Die Systemachse Z verlauft parallel zu der z-Achse 
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eines in Fig. 1 dargestellten raumlichen rechtwinkligen 
Koordinatensystems . 

Die Spalten des Detektorsystems 5 verlaufen ebenfalls in 
5 Richtung der z-Achse, wahrend die Zeilen, deren Breite b in 
Richtung der z-Achse gemessen wird und beispielsweise 1 mm 
betragt, quer zu der Systemachse Z bzw. der z-Achse 
verlaufen . 

10 Urn ein Untersuchungsob j ekt , z.B. einen Patienten, in den 

Strahlengang des Rontgens trafilenbiindels bringen zu konnen, 
ist eine Lagerungsvorrichtung 9 vorgesehen, die parallel zu 
der Systemachse Z, also in Richtung der z-Achse verschiebbar 
ist, und zwar derart, dass eine Synchronisation zwischen der 

15 Rotationsbewegung des Drehrahmens 7 und der 

Trans 1 a t i onsbewegung der Lagerungsvorrichtung in dem Sinne 
vorliegt, dass das Verhaltnis von Translations- zu 
Rotationsgeschwindigkeit konstant ist, wobei dieses 
Verhaltnis einstellbar ist, indem ein gewunschter Wert fur 

2 0 den Vorschub v der Lagerungsvorrichtung pro Umdrehung des 
Drehrahmens gewahlt wird. 

Es kann also ein Volumen eines auf der Lagerungsvorrichtung 9 
befindlichen Untersuchungsob j ekts im Zuge einer 

2 5 Volumenabtastung untersucht werden, wobei die 

Volumenabtastung in Form einer Spiralabtastung in dem Sinne 
vorgenoiranen werden kann, dass unter gleichzei tiger Rotation 
der Messeinheit 1 und Translation der Lagerungsvorrichtung 9 
mittels der Messeinheit pro Umlauf der Messeinheit 1 eine 

3 0 Vielzahl von Projektionen aus verschiedenen 

Projektionsrichtungen aufgenommen wird. Bei der 
Spiralabtastung bewegt sich der Fokus F der 

Rontgenstrahlenquelle relativ zu der Lagerungsvorrichtung 9 
auf einer in Fig. 1 mit S bezeichneten Spiralbahn. Die 
3 5 Spiralabtastung muss sich in a-Richtung iiber wenigstens 

T+2/^ax erstrecken, urn die vollstandige Rekons truktion eines 
CT-Bildes pro Zeile von Detektorelementen zu gestatten, sie 
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kann innerhalb der technischen Grenzen des CT-Gerats aber 
auch beliebig langer sein. 

Ein Volumen des Untersuchungsob j ekts kann jedoch infolge des 
Umstandes, dass mehrere Zeilen von Detektorelementen 4 
vorhanden sind, auch im Zuge einer sogenannten 
Tomogrammabtastung untersucht werden, bei der keine 
Relativbewegung in Richtung der z-Achse zwischen Messeinheit 
1 und Lagerungsvorrichtung 9 stattfindet (v=0). Im Falle der 
Tomogrammabtastung ist also die GroSe des untersuchten 
Volumens durch die Zahl der aktiven Zeilen von 
Detektorelementen 4 bestimmt. Wahrend einer 
Tomogrammabtastung bewegt sich der Fokus F auf einer 
kreisformigen Fokusbahn, die in einer im Folgenden als 
Mittelebene bezeichneten Ebene liegt. 

Die Tomogrammabtastung kann in Form eines Teilumlaufs oder in 
Form eines Vollumlaufs erfolgen, wobei der Teilumlauf ein 
Teilumlauf intervall von wenigstens /r+2#nax (ein Halbumlauf 
plus Facherof f nungswinkel ) umfasst, das eine vollstandige 
Rekonstruktion eines CT-Bildes gestattet, wahrend ein 
Vollumlauf 2zr umfasst. 

Die wahrend der Spiral- oder Tomogrammabtastung aus den 
Detektorelementen jeder aktiven Zeile des Detektorsys terns 5 
parallel ausgelesenen, den einzelnen Projektionen P(a,/3,q) in 
Facherstrahlgeometrie entsprechenden Messdaten werden in 
einer Datenauf berei tungseinheit 10 einer Digital /Analog-Wand- 
lung unterzogen, serialisiert und an einen Bildrechner 11 
ubertragen. 

Nach einer Vorverarbeitung der Messdaten in einer 
Vorverarbeitungseinheit 12 des Bildrechners 11 gelangt der 
resultierende Datenstrom zu einer 

Schnittbildrekonstruktionseinheit 13, die aus den Messdaten 
Schnittbilder von gewunschten Schichten des 
Untersuchungsob j ekts nach einem noch im Einzelnen zu 
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beschreibenden erf indungsgemaSen Verfahren auf Basis der 
"Filtered Backpro jection" rekonstruiert. 

Die CT-Bilder setzen sich aus matrixartig zusammengesetzten 
Pixeln (Pixel= picture element) zusammen, wobei die Pixel der 
jeweiligen Bildebene zugeordnet sind, jedem Pixel eine CT- 
Zahl in Hounsfield Units (HU) zugeordnet ist und die 
einzelnen Pixel entsprechend einer CT-Zahl/Grauwertskala in 
einem ihrer jeweiligen CT-Zahl entsprechenden Grauwert 
dargestellt werden. Dabei veranschaulicht jedes Pixel ein 
Voxel (Voxel = volume element) der in dem CT-Bild 
veranschaulichten Schicht des Untersuchungsob j ekts . Da 
infolge der Mehrzeiligkei t des Detektorsystems 5 und 
gegebenenfalls der Spiralabtastung Messdaten bezuglich 
mehrerer Schichten des Untersuchungsob j ekts gewonnen werden, 
stehen 3D-Daten zur Verfugung, die im Rahmen der Erfindung 
einer 3D-Ruckpro jektion unterzogen werden. Als Endergebnis 
stehen 3D-Bilddaten in Form einer dreidimensionalen Matrix, 
beispielsweise mit den Achsen x, y, z zur Verfugung, wobei 
jedes Element der Matrix einem Voxel V(x,y,z) entspricht und 
den der zugehorigen CT-Zahl entsprechenden Grauwert enthalt. 
Diejenigen Elemente der dreidimensionalen Matrix die den 
gleichen x- , y- , oder z-Wert aufweiseh, stellen dann jeweils 
ein planares Schnittbild der dem ma&geblichen x- , y- , oder z- 
Wert entsprechenden Schicht des Untersuchungsob j ekts dar . 

Die von der Schnittbildrekonstruktionseinhei t 13 
rekonstruierten Bilder werden auf einer an den Bildrechner 13 
angeschlossenen Anzeigeeinhei t 16, z.B. einem Monitor, 
dargestellt . 

Die Rontgenstrahlenquelle 2, beispielsweise eine 
Rontgenrohre, wird von einer Genera tor einheit 17 mit den 
notwendigen Spannungen und Stromen, beispielsweise der 
Rohrenspannung U, versorgt. Urn diese auf die jeweils 
notwendigen Werte einstellen zu konnen, ist der 
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Generatoreinheit 17 eine Steuereinheit 18 mit Tastatur 19 
zugeordnet, die die notwendigen Eins tellungen gestattet. 

Auch die sonstige Bedienung und Steuerung des CT-Gerates 
erfolgt mittels der Steuereinheit 18 und der Tastatur 19, was 
dadurch veranschaulicht ist, dass die Steuereinheit 18 mit 
dem Bildrechner 11 verbunden ist. 

Unter anderem kann die Anzahl Q der aktiven Zeilen von 
Detektorelementen 4 und damit die Position der 
Strahlenblenden 3 und 6 eingestellt werden, wozu die 
Steuereinheit 18 mit den Strahlenblenden 3 und 6 zugeordneten 
Verstelleinheiten 20 und 21 verbunden ist. Weiter kann die 
Rotations zeit r eingestellt werden, die der Drehrahmen 7 fur 
eine vollstandige Umdrehung benotigt, was dadurch 
veranschaulicht ist, dass die dem Drehrahmen 7 zugeordnete 
Antriebseinheit 22 mit der Steuereinheit 18 verbunden ist. 

Obwohl es grundsatzlich moglich ist, das erf indungsgemaSe 
Verfahren auch in Fachers trahlgeometrie zu realisieren, wird 
das beschriebene CT-Gerat vorzugsweise in einem Modus 
betrieben, in dem das erf indungsgemaSe Verfahren in 
Parallelstrahlgeometrie realisiert ist. 

Demnach werden die bei der Abtastung des fur die jeweilige 
Untersuchung relevanten Korperbereichs des Patienten durch 
Spiral- Oder Tomogrammabtas tung in Facherstrahlgeometrie 
gewonnenen Daten zunachst in an sich bekannter Weise durch 
ein im Allgemeinen als "rebinning" bezeichnetes Verfahren in 
Daten in Parallelstrahlgeometrie umgewandelt. Diese 
Umwandlung beruht auf einer Umsortierung der in 
Facherstrahlgeometrie gewonnenen Daten derart, dass aus 
unterschiedlichen in Facherstrahlgeometrie auf genommenen 
Projektionen Strahlen entnommen und zu einer Projektion in 
Parallelstrahlgeometrie zusammengef ugt werden. In 
Parallelstrahlgeometrie reichen Daten aus einem Intervall der 
Lange n aus, urn ein volls tandiges Bild rekons truieren zu 
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konnen. Um diese Daten gewinnen zu konnen, rtriissen 
nichtsdestotrotz zur Daten in Facherstrahlgeometrie aus einera 
Intervall der Lange ^+2y5 ma x zur Verfugung stehen. 

Eine Projektion in Parallelstrahlgeometrie ist in Fig. 3 
veranschaulicht. Demnach nehmen alle n Parallelstrahlen RPi 
bis RP N dieser Projektion gegenuber der x-Achse des in Fig. 3 
darges tell ten und mit dem gemaE Fig. 1 iibereinstimmenden 
Koordinatensystem den Parallelf acherwinkel 0 ein. 

Anhand des in Fig. 3 in ausgezogener Linie dargestellten 
Parallelstrahls RPi soil im Folgenden der Ubergang von 
Facherstrahl- auf Parallelstrahlgeometrie erlautert werden. 

Der Parallelstrahl RPi entstammt einer fur die auf der 
Fokusbahn S liegende Fokusposition Fi in 

Facherstrahlgeometrie gewonnenen Projektion. Der zu dieser 
Projektion in Facherstrahlgeometrie gehorige, durch die 
Drehachse 14 und damit die z-Achse des Koordinatensys terns 
verlaufende Zentralstrahl RF Z1 ist in Fig. 3 ebenfalls 
eingetragen. Der . Fokusposition Fi entspricht der Fokuswinkel 
an dies ist der Winkel , den die x-Achse und der 
Zentralstrahl RF Z i einschliefcen . Gegenuber dem Zentralstrahl 
RFzi weist der Strahl RPi den Facherwinkel (5 auf. Somit ist 
leicht erkennbar, dass fur den Parallelf acherwinkel 6 gilt: 
0 = a + J3. 

Der rechtwinklig zum jeweiligen Parallelstrahl gemessene 
Strahlabstand p von der Drehachse 14 bzw. der z-Achse ist 
durch p=R F sin (fi) gegeben . 

Wie anhand des in Fig. 3 in verstarkter Linie dargestellten, 
durch die Drehachse 14 bzw. die x-Achse verlaufenden 
Zentralstrahls RP Z deutlich wird, handelt es sich bei diesem 
Strahl um den Zentralstrahl einer in Fachergeometrie fur die 
Fokusposition F z unter dem Fokuswinkel a z aufgenommene 
Projektion in Facherstrahlgeometrie. Da fur den Zentralstrahl 
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einer in Facherstrahlgeometrie auf genommenen Projektion J3=0 
gilt, wird deutlich, dass fur den Fall von Zentrals trahlen 
gilt : 

5 Je nachdem, ob ein azimutales oder volls tandiges "rebinning" 
durchgefuhrt wird, liegen die Parallelpro j ektionen in der 
Form 

P(a,P,q) 
oder in der Form 
10 P{0,P,q) 
vor, wobei 

a der Fokuswinkel 

P der Facherwinkel 

q der der z-Koordinate entsprechende Zeilenindex 

15 des Detektorsys terns , 

Q = a + P der Parallelf acherwinkel 

p = R F sin(p) die dem Strahlabs tand von der Drehachse 

(Systemachse) entsprechende Parallelkoordinate , 
und 

2 0 Rf der Radius der Fokusbahn sind. 

In einer ersten Ausf uhrungsf orm des erf indungsgemafien 
Verfahrens entsprechenden, mittels der Tastatur 19 wahlbaren 
ersten Betriebsart arbeitet das beschriebene CT-Gerat auf 

2 5 Basis von durch azimutales "rebinning" gewonnenen 

Pro j ektionen. Die diesen Pro j ektionen entsprechenden Daten 
werden im Falle des azimutalen "rebinnings" in /?-Richtung 
gefiltert, d.h. jeweils in Richtung der zu der Fokusposi tion 
des Zentralstrahls der jeweiligen Parallelpro j ektion 

30 gehorigen Tangente T (siehe Fig. 3), und zwar unter 

Verwendung eines der in der Computer tomographie ublichen 
Filterkerne, z.B. Shepp-Logan - oder Ramachandran- 
Lakshminarayanan-Kern . 



35 



Die derart gefilterten Paralleldaten werden dann in der in 
der Weise ruckpro j iziert , dass im Zuge der Riickproj ektion fur 
jedes Voxel V(x,y,z) fur jedes 0 <=[0,7r[ fur die Strahlen 
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p(& + kn, fi, q) , deren Projektion entlang der Systemachse durch 
(x,y) geht, die Summe 

gebildet wird, wobei 

x.y.z die Koordinaten des jeweiligen Voxels V(x,y,z) sind, 

A: eine ganze Zahl entsprechend der Anzahl der in die 

Rekonstruktion einbezogenen Halbumlaufe des Fokus , 
P die Facherwinkel derjenigen Strahlen sind, deren 

Projektionen langs der Systemachse durch die 

Koordinaten (x,y) jeweiliges Voxels V(x,y,z) verlaufen, 

und 

h eine die Schichtdicke der in dem erzeugten 

Schnittbild dargestellten Schicht des 

Untersuchungsobjekts bestimmende Gewichtungsf unktion 
ist, sowie 

d eine Funktion ist, die gleich dem Abstand des 

jeweiligen Strahls von dem entsprechenden Voxel V(x,y) 
ist bzw. von dem Abstand des jeweiligen Strahls von 
dem entsprechenden Voxel V(x,y) abhangt, und 

W(q) eine Wichtungs funktion darstellt, die Strahlen, die 
im Strahlenbundel bezogen auf die Ausdehnung in 
Rotationsachsenrichtung peripher oder zumindest 
dezentral positioniert sind, geringer wichtet als 
Strahlen, die zentral angeordnet sind. 

infolge der gewahlten Filterrichtung und der infolge der 
Summation sowohl uber k als auch uber q werden zum einen cone 
beam artifacts vermieden, im Interesse einer hohen 
Dosisnutzung alle durch ein Voxel V(x,y,z) verlaufende Strahlen 
berucksichtigt, jedoch zusatzlich auch die Posit ionierung des 
Strahls im Strahlenbundel durch entsprechende Wichtung 
berucksichtigt . 
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Der einem Voxel V(x,y,z) zugeordnete Absorptionswert /u xyz wird 
durch Summation uber 6 uber mindestens einen Halbumlauf 
gewonnen, d.h. durch Bildung von 

^,,.*=Z^» 

e 

Die dem jeweiligen Absorptionswert entsprechende CT-Zahl wird 
in herkommlicher Weise aus dem Absorptionswert ermittelt. 

Dabei konnen mittels der Tastatur 19 unterschiedliche 
Gewichtungsfunktionen h und unterschiedliche Funktionen d 
eingestellt werden. 

Als Gewichtungsfunktion h eignet sich beispielsweise eine 
Dreiecks- oder Trapezf unktion . 

Die Gewichtungsfunktion W ist als eine Funktion des 
Projektionswinkels 6 dargestellt. 

Als Funktion d kann der Abstand des jeweiligen 
Parallelstrahls vom Voxel V(x t y,z) oder stattdessen 
beispielsweise die z(axiale) -Komponente dieses Abstandes 
eingestellt werden . 

In einer Abwandlung der vorstehend beschriebenen ersten 
Betriebsart wird zur Ruckpro j ektion der Paralleldaten die auf 
die Summe H der Gewichte h 

h = j Z ^fe) • h ( d *^ ( 0 + k *> *)) 



normierte Summe 



± Z Z w b) ■ h ( d *-y* b + k "> & *))• p ( 0 +k *> /3 ' q } 

H k q 
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gebildet. Dies ermoglicht eine nochmals verbesserte 
Bildqualitat , da eine mogliche Uberbetonung von Voxeln, die 
in mehreren Halbumlaufen "beleuchtet " , also von Strahlen 
getroffen werden, beseitigt wird und somit entsprechende 
Artefakte vermieden werden. Diese Redundanz tritt bei 
Spiralabtastungen dann auf , wenn die pro Vollumlauf der 
Messanordnung stattf indende Relatiwerschiebung so gering ist 
(geringer Pitch), dass Voxel mehrfach bestrahlt werden. 

Eine einer weiteren Ausf uhrungsf orm des erf indungsgemaSen 
Verfahrens entsprechenden, mittels der Tastatur 19 wahlbaren 
zweite Betriebsart unterscheidet sich von der ersten 
Betriebsart dadurch, dass das beschriebene CT-Gerat nicht auf 
Basis von durch azimutales, sondern von durch vollstandiges 
"rebinning" gewonnenen Projektionen arbeitet. Die diesen 
Projektionen entsprechenden Daten werden im Falle des 
vollstandigen "rebinnings" in p-Richtung gefiltert, also 
jeweils ebenfalls in Richtung der zu der Fokusposition des 
Zentralstrahls der jeweiligen Parallelpro j ektion gehorigen 
Tangente T (siehe Fig. 3) . 

Demnach wird fur die derart gefilterten Paralleldaten im Zuge 
der Ruckproj ektion die Summe 

k q 

gebildet, wobei p die Parallelkoordinaten derjenigen 
Strahlen sind, deren Projektionen langs der Systemachse durch 
die Koordinaten fay) jeweiliges Voxels Vfay.z) verlaufen. 

Auch im Falle der zweiten Betriebsart wird in einer 
Abwandlung zur Ruckproj ektion der Paralleldaten eine auf die 
Summe H der Gewichte h 

H = £ £ W(q) ■ h{d w (0 + kx,p,q)) 
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normierte Summe, namlich die Summe 

Ky, (^) = 77 Z Z * ^ + ^ + ^ * ^ 

H k q 

gebildet . 

Im Falle der vorstehend beschriebenen ersten und zweiten 
Betriebsarten wird eine Funktionsweise des erf indungsgemaSen 
Verfahrens vorgesehen, bei dem im Zusammenhang mit einem 
Voxel V(x,y,z) alle Strahlen betrachtet werden, deren 
Projektion entlang der Drehachse 14 bzw. der z-Achse durch x, 
y geht. Ob und in welchem MaSe diese Strahlen berucksichtigt 
werden, wird durch die Gewichtungsf unktionen W und h, sowie 
die Funktion d bestinunt. 

Das CT-Gerat kann aber auch mittels der Tastatur 19 wahlbare 
weitere Betriebsarten aufweisen, die den zuvor beschriebenen 
mit dem Unterschied entsprechen, dass fur eine gegebene 
Fokusposition der durch das jeweilige Voxel V(x,y,z) 
verlaufende theoretische Strahl ermittelt wird und dann unter 
Berucksichtigung der Gewichtungsf unkt ion h und der Funktion d 
in die Summenbildung im Zuge der Ruckpro j ektion nur 
diejenigen Strahlen einbezogen werden, die tatsachlich einen 
Beitrag zu der Summe liefern konnen. 

Im Falle der beschriebenen Ausf iihrungsbeispiele wird die 
Relativbewegung zwischen der Messeinheit 1 und 
Lagerungsvorrichtung 9 jeweils dadurch erzeugt, dass die 
Lagerungsvorrichtung 9 verschoben wird. Es besteht im Rahmen 
der Erfindung jedoch auch die Moglichkeit, die 
Lagerungsvorrichtung 9 ortsfest zu lassen und statt dessen 
die Messeinheit 1 zu verschieben. AuEerdem besteht im Rahmen 
der Erfindung die Moglichkeit, die notwendige Relativbewegung 
durch Verschiebung sowohl der Messeinheit 1 als auch der 
Lagerungsvorrichtung 9 zu erzeugen. 
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In der Figur 4 wird, vereinf achend und schematisch 
dargestellt, das Problem der Wichtung der Strahlen in einer 
2D-Darstellung gezeigt . In diesem Beispiel ist das Gantry des 
CT's mit Fokus und Detektor in 0°-Stellung durchgezogen und 
5 urn 180° urn die z-Achse gedreht in 180°-Stellung gepunktet 
gezeichnet. Der Detektor weist 12 Detektorzeilen LI bis L12 
auf, die jeweils von den zugeordneten Strahlen SI bis S12 des 
Strahl enbiindels getroffen werden, Es ist darauf hinzuweisen, 
dass die gezeigten raumlichen Dimensionen bezuglich der 
10 Detektorausdehnung ubertrieben dargestellt sind, urn die 
Erfindung anschaulicher erklaren zu konnen. Das gezeigte 
Gantry bewegt sich hier auf einer Kreisbahn urn die z-Achse 
oder Systemachse. Eine Bewegung in Richtung z-Achse findet 
nicht statt. 

15 

Betrachtet man das Voxel V2 , das zentral im Strahl enbiindel 
liegt, so wird dieses Voxel in beiden Stellungen der Gantry 
von den zentral liegenden Strahlen S4 durchdrungen . 
Erf indungsgemaS werden die Messungen dieser Strahlen hoch 

20 gewichtet. Gleichzeitig wird in der 0°-Stellung das Voxel VI, 
das etwas peripherer liegt, von dem auSen liegenden Strahl SI 
in der 0°-Position und dem weiter innen liegenden Strahl S4 
durchdrungen. Entsprechend der erf indungsgemaSen Wichtung 
wird dem Strahl Si weniger Gewicht beigemessen als dem Strahl 

25 S4. Betrachtet man das ebenfalls peripher angeordnete Voxel 
V3, so ist hier die Wichtung umgekehrt . Da hier der Detektor 
mit Detektorelementen aufgebaut ist, die matrixartig in 
Zeilen und Spalten aufgeteilt ist, kann die Wichtung der 
Strahlen direkt an die Zeilennummer g gekoppelt werden. 

30 

In der darunter liegenden Figur 5 sind drei Beispiele eines 
erf indungsgemaSen Verlaufs von Wichtungen W q dargestellt. 
Allen gemeinsam ist die starkere Wichtung der zentralen 
Strahlen gegeniiber den auSeren Strahlen. Die Funktionen sind 
35 allerdings unterschiedlich . Die gepunktete Linie stellt einen 
beidseits vom Zentrum linear abfallenden Verlauf dar, die 
gestrichelte Linie ergibt eine gleiche Wichtung uber einen 
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groSen zentralen Bereich mit peripher schnell abfallender 
Tendenz, und die durchgezogene Linie zeigt einen 
glockenf ormigen Verlauf , der einen weichen Ubergang zwischen 
den Extremwerten 0 und 1 erlaubt. 

5 

Die Figur 6 zeigt die Strahl enbundel einer Gantry eines 
Spiral-CT's jeweils in 0°-Position und benachbarter 180°- 
Position. Entsprechend dem Wesen des Spiral-CT's sind diese 
Positionen zusatzlich zur Drehung in z-Richtung zueinander 

10 versetzt. In dieser Darstellung wird das Voxel VI durch den 
Strahl SI sowohl in der 0°-Position als auch der 180°- 
Position durchdrungen . Entsprechend dieser dezentralen Lage 
des Strahls SI beziehungsweise der zugehorigen auSen im 
Detektor liegenden Detektorzeile Ll werden diese beiden 

15 Messungen gering bewertet . Das Voxel V2 wird hingegen in der 
0°-Position vom dezentral liegenden Strahl S12 durchdrungen, 
der in dieser Position fur das Voxel V2 gering gewichtet 
wird, in der 180°-Position jedoch vom zentralen Strahl S6 
durchdrungen wird, der fur dieses Voxel V2 hoch gewichtet 

2 0 wird. Entsprechend der Lage der Strahlen SI in der 0°- und 

180°-Position werden diese Strahlen in Bezug auf das Voxel VI 
beide gering bewertet. Insgesamt ergibt sich durch diese hier 
dargestellte zusatzliche Gewichtung W eine wesentliche 
Verbesserung der Bildqualitat . 

25 

Es ist darauf hinzuweisen, dass das konusf ormige 
Rontgenstrahlenbundel im Falle des beschriebenen 
Ausf uhrungsbeispiels einen rechteckigen Querschnitt aufweist, 
im Rahmen der Erfindung jedoch auch andere 

3 0 Querschnittsgeometrieri moglich sind. 

Im Zusammenhang mit den vorstehend beschriebenen 
Ausfiihrungsbeispielen finden CT-Gerate der dritten Generation 
Verwendung, d.h. die Rontgens trahlenquelle und das 
3 5 Detektorsystem werden wahrend der Bilderzeugung gemeinsam urn 
die Systemachse verlagert. Die Erfindung kann aber auch im 
Zusammenhang mit CT-Geraten der vierten Generation, bei denen 
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nur die Ront gens t r ahl enquel 1 e um die Systemachse verlagert 
wird und mit einem f ests tehenden Detektorring zusammenwirkt, 
Verwendung finden, sofern es sich bei dem Detektorsys tern um 
ein mehrzeiliges Array von Detektorelementen handelt. 

5 

Auch bei CT-Geraten der funften Generation, d.h. CT-Geraten, 
bei denen die Rontgenstrahlung nicht nur von einem Fokus , 
sondern von mehreren Foken einer oder mehrerer um die 
Systemachse verlagerter Rontgens trahlenquellen ausgeht, kann 
10 das erfindungsgemaSe Verfahren Verwendung finden, sofern das 
Detektorsystem ein mehrzeiliges Array von Detektorelementen 
aufweist . 

Die im Zusammenhang mit den vorstehend beschriebenen 
15 Ausf uhrungsbeispielen verwendeten CT-Gerate weisen ein 
Detektorsystem nach Art einer orthogonalen Matrix 
angeordneten Detektorelementen auf . Die Erfindung kann aber 
auch im Zusammenhang mit CT-Geraten Verwendung finden, der en 
Detektorsystem in einer anderen Weise ein flachenhaft 
20 angeordnetes Array mit Detektorelementen aufweist. 

Die vorstehend beschriebenen Ausf iAhrungsbeispiele betreffen 
die medizinische Anwendung des erfindungsgemaSen Verf ahrens . 
Die Erfindung kann jedoch auch aufierhalb der Medizin, 
25 beispielsweise bei der Gepackpruf ung oder bei der 
Ma terialunters uc hung , Anwendung finden. 



200210627 



23 

Patent anspruche 

1. Verfahren zur Bilderstellung in der 
Computertomographie, aufweisend die Verf ahrensschri tte : 

5 1.1. zur Abtastung eines Untersuchungsobjekts mit einem von 
mindestens einem Fokus ausgehenden Strahlenbundel und 
mit einem flachig ausgebildeten Detektorarray mit einer 
Vielzahl von verteilten Detektorelementen zum 
Detektieren der Strahlen des Strahlenbiindels wird der 

10 mindestens eine Fokus relativ zu dem 

Untersuchungsobjekt auf mindestens einer das 
Untersuchungsobjekt umlaufenden Fokusbahn mit 
gegeniiberliegendem Detektorarray rotiert, wobei die 
Detektorelemente des Detektorarray Ausgangsdaten 

15 liefern, welche die Schwachung der Strahlen beim 

Durchgang durch das Untersuchungsobjekt reprasentieren, 

1.2. die Ausgangsdaten werden gefiltert, 

1.3. die gefilterten Ausgangsdaten werden anschlieSend zur 
Erzeugung wenigstens eines Schnittbildes einer eine 

2 0 Schichtdicke aufweisenden Schicht des 

Untersuchungsobj ekts dreidimensional riickpro j iziert , 
wobei das Schnittbild aus den Ausgangsdaten gewonnene 
Absorptionswerte der zu der Schicht des 

Untersuchungsobjekts gehorigen Voxel fur die Strahlung 
2 5 des Strahlenbiindels reprasentiert, und 

1.4. bei der Riickpro j ektion die Strahlen in Abhangigkeit von 
ihrer Position im Strahlenbundel gewichtet werden. 

2. Verfahren nach Anspruch 1, dadurch 

30 gekennzeichnet, dass die Filterung in Richtung 

der zu der jeweiligen Fokusposi tion gehorigen Tangente 
an die Fokusbahn vorgenommen wird. 

3. Verfahren nach einem der Anspruche 1 bis 2, dadurch 
35 gekennzeichnet, dass das Strahlenbundel eine 

Ausdehnung in Rotationsrichtung und eine Ausdehnung in 
Rotationsachsenrichtung aufweist, wobei Strahlen, die 
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im Strahlenbundel bezogen auf die Ausdehnung des 
Strahlenbundels in Rotationsachsenrichtung zentral 
angeordnet sind, starker gewichtet werden, als die 
Strahlen, die im Strahlenbiindel randnah bezogen auf die 
5 Ausdehnung des Strahlenbundels in 

Rotationsachsenrichtung angeordnet sind. 

4 . Verf ahren nach einem der Anspruche 1 bis 3, dadurch 
gekennzeichnet , dass vor dem Filtern eine 
10 Umrechnung der in Facherstrahlgeometrie gewonnenen 

Ausgangsdaten in Form von Strahlen P(a,/3,q) in in 

Parallelstrahl-Geometrie vorliegende Paralleldaten in 
Form von Strahlen P(&,j3,q) bzw. P(0,p r q) vorgenommen 

wird, wobei 
15 a der Fokuswinkel , 

P der Facherwinkel , 

q der der z-Koordinate entsprechende 

Zeilenindex des Detektorsys terns , 
0 = a + j3 der Parallel f acherwinkel , 

20 p=R F sin(p) die dem Strahlabstand von der Drehachse 

(Systemachse) entsprechende 
Parallelkoordinate , und 
R F der Radius der Fokusbahn sind. 



25 



30 



Verf ahren nach Anspruch 4, dadurch 
gekennzeichnet, dass die Riickpro jektion der 
Paralleldaten in der Weise erfolgt, dass im Zuge der 
Ruckpro j ektion fur jedes Voxel V(x,y t z) fur jedes #e[0,7r[ 

fur die Strahlen p(& + kn, fi,q) bzw. P{p + kn,p^q) , der en 
Pro jektion entlang der Systemachse durch (x,y) geht, die 
Summe 



x r y 7 z 



0 + kx,^,q 



gebildet wird, wobei 

x.y.z die Koordinaten des jeweiligen Voxels V(x,y,z) sind, 
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k eine ganze Zahl entsprechend der Anzahl der in 

die Rekonstruktion einbezogenen Halbumlaufe des 
Fokus , 

p die Parallelkoordinaten derjenigen Strahlen sind, 

5 deren Projektionen langs der Systemachse durch 

die Koordinaten (x,y) jeweiliges Voxels V(x,y,z) 
verlauf en, 

P die Facherwinkel derjenigen Strahlen sind, deren 

Projektionen langs der Systemachse durch die 
10 Koordinaten (x,y) jeweiliges Voxels V(x,y,z) 

verlauf en, und 
h eine die Schichtdicke der in dem erzeugten 

Schnittbild darges tell ten Schicht des 
Untersuchungsob j ekts bes timmende 
15 Gewichtungsf unktion ist, sowie 

d eine Funktion ist, die gleich dem Abstand des je- 

weiligen Strahls von dem entsprechenden Voxel 
V(x,y) ist bzw. von dem Abstand des jeweiligen 
Strahls von dem entsprechenden Voxel V(x,y) 
2 0 abhangt , und 

W eine Wichtungs funktion darstellt, welche Strahlen 
mit einem betragsmaSig groSen 
Parallelf acherwinkel 6 geringer wichtet als 
Strahlen mit einem kleinem Parallelf acherwinkel 
25 0. 



30 



Verfahren nach Anspruch 5, dadurch 

gekennzeichnet, dass bei der Ruckpro j ektion der 
Paralleldaten die auf die Sunune H der Gewichte h 



x, y.z 



0 + knA ~\,q 



normierte Sunune 

H k q 



gebildet wird. 
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7. Verfahren nach einem der Anspruche 5 bis 6, dadurch 
gekennzeichnet, dass die Wichtungsf unktion eine 
Funktion des Parallelf acherwinkels mit W ( 9 + kn ) 
darstellt . 

5 

8. Verfahren nach Anspruch 7 , dadurch 
gekennzeichnet , dass die Wichtungsf unktion 
VJ(0 + kn) eine glatte Funktion darstellt, die fur 
zentral angeordnete Strahlen im Strahlenbundel bezogen 

10 auf den Parallelf acherwinkel den Wert 1 aufweist und 

fur randseitig angeordnete Strahlen gegen Null geht . 

9. Verfahren nach einem der Anspruche 4 bis 5, dadurch 
gekennzeichnet , dass das Detektorarray 

15 zeilenartig angeordnete Detektorelemente aufweist und 

die Wichtungsf unktion eine Funktion der Zeilennummer 
W(q) darstellt. 



10. Verfahren nach Anspruch 8, dadurch 

20 gekennzeichnet, dass die Wichtungsf unktion W(g) 

eine glatte Funktion darstellt, die fur Strahlen zu 
mindestens einer zentral liegenden Detektorzeile den 
Wert 1 aufweist und fur Strahlen zu randnahen 
Detektorzeilen gegen Null geht. 

25 

11. Verfahren nach einem der Anspruche 1 bis 9, dadurch 
gekennzeichnet , dass die Fokusbahn eine 
Kreisbahn ist. 



3 0 12. Verfahren nach einem der Anspruche 1 bis 9, dadurch 
gekennzeichnet , dass die Fokusbahn eine 
Spiralbahn ist, die dadurch zustande kommt, dass der 
Fokus auf einer Kreisbahn urn die Systemachse bewegt 
wird und gleichzeitig eine Relativbewegung zwischen 

3 5 Fokus und Untersuchungsob j ekt in Richtung der 

Systemachse erf olgt . 
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13. Verfahren nach einem der Anspruche 1 bis 12, dadurch 
gekennzeichnet , dass die Detektorelemente auf 
dem Detektorarray in Zeilen und Reihen verteilt, 
vorzugsweise matrixartig, angeordnet sind. 

14. CT-Gerat zur Abtastung eines Untersuchungsobjekts mit 
einem von mindestens einem Fokus ausgehenden 
Strahlenbiindel und mit einem flachig ausgebildeten 
Detektorarray mit einer Vielzahl von verteilten 
Detektorelementen zum Detektieren der Strahlen des 
Strahlenbiindels , wobei der mindestens eine Fokus 
relativ zu dem Untersuchungsob j ekt auf mindestens einer 
das Untersuchungsob j ekt umlaufenden Fokusbahn mit 
gegeniiberliegendem Detektorarray bewegt wird, wobei 
zumindest Mittel zur Sammlung von Detektordaten, 
Filterung und Riickpro j ekt ion gemaS einem der Anspruche 
1 bis 13 vorgesehen sind. 

15. CT-Gerat gemafi Anspruch 14, dadurch 
gekennzeichnet , dass die genannten f unktionalen 
Mittel zumindest teilweise durch Programme oder 
Programm-Module verwirklicht werden . 
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Zusammenf assung 

Verfahren zur Bilderstellung in der Computertomographie und 
CT-Gerat zur Durchfiihrung des Verfahrens 

5 

Die Erfindung betrifft ein Verfahren zur Bildrekonstruktion 
fur die Computertomographie mit einem flachig ausgedehnten 
Detektor, wobei bei der Riickpro j ektion die Strahlen in 
Abhangigkeit von ihrer Position im Strahlenbundel gewichtet 
10 werden. 

AuSerdem betrifft die Erfindung ein CT-Gerat zur Durchfuhrung 
des oben genannten Verfahrens. 

15 FIG 4 
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Bezugszeichenliste 



1 Messanordnung / Messeinheit 

2 Rontgens trahlenquel 1 e 

5 3 quellennahe Strahlenblende 

4 Detektorelemente 

5 Detektorsystem 

6 detektornahe Strahlenblende 

7 Drehrahmen 
10 8 Rands trahl en 

9 Lagervorr i chtung 

10 Datenauf bereitungseinhei t 

11 Bildrechner 

1 2 Vorverarbei tungseinhei t 

15 13 Schni t tbi ldr ekons trukt i ons einhei t 

14 Drehachse 

16 Anzeigeeinheit 

17 Genera toreinheit 

18 Steuereinheit 

20 19 Tastatur der Steuereinheit 

2 0 Verstelleinheit 

21 Verstelleinheit 

22 Antriebseinrichtung 
Sx Strahlen 

25 Lx Detektorzeilen 

Vx Voxel 

F Fokus 

F x Fokusposition 

S Spiralbahn 
3 0 Z Systemachse 

U Rohrenspannung 

a Fokuswinkel 

J3 Facherwinkel 

q Zeilenindex der z-Koordinate 

35 0 = a + /3 Parallel facherwinkel 



p = R F sin(j3) Parallelkoordinate des Strahlanstandes der 
Drehachse 



Radius der Fokusbahn 

Koordinaten des jeweiligen Voxels V(x,y,z) 
ganze Zahl entsprechend der Anzahl der in die 
Rekonstruktion einbezogenen Halbumlaufe des Fokus 
Parallelkoordinaten derjenigen Strahlen, deren 

Projektionen langs der Systemachse durch die 
Koordinaten (x t y) des jeweiligen Voxels (x t y,z) 
verlauf en 

Facherwinkel derjenigen Strahlen, deren 
Projektionen langs der Systemachse durch die 
Koordinaten (x,y) des jeweiligen Voxels (x,y,z) 
verlauf en 

schichtdickenbestimmende Gewichtungsf unktion 
Funktion der Abstandes des jeweiligen Strahls vom 
Voxel 

Wichtungs funktion fur die Position des Strahls im 
S t r ah 1 enbiinde 1 



